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The performance of a silicon device/l.C. is dominantly 
influenced by the impurity distribution in the active regions 
of the device. These impurities are introduced into silicon 
by diffusion/ion implantation. In the device/l.C. fabrication 
process sequence, this impurity deposition is invariably 
followed by other high temperature processes such as drive-in, 
oxidation etc. The deposited impurities redistribute themselvs 
during these processes and the final impurity distri-^ 
bution is rather a complicated function of the initial distri- 
bution and the parameters such as ambient, temperature etc., oi 
the later processes. To accurately relate the process para- 
meters with the device performance, therefore, a precise know- 
ledge of the effect of process parameters on the impurity diffu 

j 

sion is essential. The recent years have soon considerable 
activity in this area, primarily motivated by two factors. 

(i) As the I.C. complexity increases, the computer aided desigil 

methods are to be used, in order to reduce the time required 
for the development of a new product. This in turn needs 


accurate mathematical models for process variable-impurity 
profile relationship, besides the models for dovices/circuits. 
(ii) It is recognised that the process-profile rolationship is 
still not adequately known and remains a major limiting factor 
in our capability to accurately relate the process variables 
with device/circuit performance. 

The objective of this work is to study tho effect of 
oxygen and chlorine (in diffusion ambient) on the diffusion of 
phosphorus in silicon. It is now well established that the 
thermal oxidation of silicon causes an injection of excess self- 
interstitials into tho bulk silicon, which cause an enhancement 
in impurity diffusion coefficient. Tho study of this oxidation 
enhanced diffusion (OED) has boon reported by several workers 
[1-3] f who estimated OED by evaluating the redistribution of 
impurities, deposited by ion implantion/CVD, during oxidation* 
Tho extraction of the diffusion coefficient in these experiments 
is todius and tho finer details of tho experimental profile 
vis-a-vis its theoretical counterpart are likely to bo ovorlooke 
duo to the approximations made for initial distribution* Tho 
effect of chlorine on phosphorus diffusion has boon reported in ; 
detail only by Nabota et al. [4]. Their experimental results, 
however, cover only a narrow range of tonporature and have boon 

OA i 

questioned duo to the high surface concentrations (> 10 cm ) 
involved in the experiment [5]. In tho present work, the effect 
of ambient has been studied by carrying out constant source 



diffusion using PHg in different ambients, under intrinsic 
diffusion conditions. This exporirriuntal situation is amenable 
to simple mathematical analysis. In addition, the largo segre- 
gation coefficient of phosphorus in Si 02 -Si system, allows the 
diffusions to bo carried out in strongly oxidising ambient, thus 
permitting the observation of OED of nearly the same magnitude 
as observed during oxidation in dry O 2 . 

A mathematical model for impurity diffusion in oxidising 
ambient from a constant source is presented. The diffusion 
equations in dopant silicate glass (DSG) and silicon are solved 
for impurity distribution in the two mediums. A reaction rate 
limited flux boundary condition, instead of that of thermo- 
dynamic equilibrium, has boon assumed at the interface. The 
general solution for impurity distribution in silicon is found ‘ 
to be a scries, which converges reasonably fast only for small 
diffusion lengths. For extreme cases such as negligible EBG 
growth or largo values of reaction rate constant/diffusion 
time, the model yields solutions, identical to those reported 
by Smits [’6] and Grove [?]• 

Phosphorus was diffused in <111> and <1CX3> silicon, using ' 
PH^ source (1000 ppmPH^^ in Ar) at various temperatures in the 
range 950°C-1250°C, In throe different ambients; (i) 95.5 W- Ng- 

1.5 W* 0-5 + 3.0 W* (PHo+Ar), (ii) 90.1 W 0^ + 9.9 W* (PH-^+Ar) 
and (iii) 94.5 W O 2 + 0*5 Vyi (PH 3 +Ar) + 5,0 Vj4 N 2 bubbled 



through TCE at room temperature. The impurity distribution 
in each sample was determined from sheet resistivity measure- 
ments made after successive removal of silicon by anodic oxida- 
tion. The diffusion parameters were then calculated by compar- 
ing the experimental profile with the appropriate theoretical 
profile derived in the model. 

In agreement with the theoretical predictions for diffusion 
in oxidising ambient, the surface impurity concentration is 
observed to be a strong function of the ratios of the PSG thick- 
ness to impurity diffusion lengths in PSG and silicon respecti- 
vely. For identical diffusion ambients, the oxidation rate is 
seen to cause the surface impurity concentration to vary by two 
orders of magnitude. On the basis of the model predictions and 
the experimental results, it is concluded that the desirable 
process features for low concentration phosphorus diffusion are 
high temperature and strongly oxidising ambient. 

In strongly oxidising ambient, the observed OED is in good 
agreement with the data reported by Antoniadis ot al. [3], for 

OED during dry O 2 oxidation. In the present work the surface 

18 

impurity concentration was of the order of 10 cm . At these 
concentrations, the OED was observed to be independent of cry- 
stal orientation and diffusion time. Various models proposed 
for the kinetics of excess-self interstitials generated during 
oxidation, are examined in the light of the observations of this 
work coupled with the observations on other oxidation related 



phenomena, such as growth of oxidation induced stacking faults 
(OSF), reported in litoraturo. It is shown that only the model 
proposed by Hu [8], satisfactorily explains all the experimental 
observations reported so far on OED and OSF growth in silicon. 

It is suggested that the surface regrowth process in Hu’s model 
depends upon the impurity typo (through its segregation property 
in SlOj-Si system) and concentration at the interface. 

The presence of chlorine in diffusion ambient is observed 
to cause a reduction in diffusion coefficient at all the tom- 
poraturcs. This effect is significantly larger on <111> than on 
<100> silicon and decreases with temperature, becoming negli- 
gible at 950°C. The reduction in diffusion coefficient is con- 
sistent with the observation of retarded growth of OSF in Cl2“02 
ambient [9], reported in literature and is a pointer to the fact 
that the presence of chlorine in the oxidation/diffusion andaient 
reduces the not generation of excess self interstitials at the 
interface. The experimental data of this work have boon corre- 
lated with the reported segregation properties of chlorine in 
Si02-Si system, during chlorine-oxidation [lO-ll] and it is 
confirmed that the chlorine containing interface acts as a sink 
for self-interstitials. Hu’s model is extended to take into 
account the effect of chlorine on diffusion coefficient. 



Tho last but important observation is the reduction o. 
diffusion coefficient in chloro-oxidising ambient below its 
value in inert ambient (by about at 1250'^C) in <111> 
silicon at high temperatures (>1150°C), It is argued to bo 
a direct ovidonco to the fact that in inert ambient , phosphorus 
diffusion proceeds via partial/pure intorstitialcy mechanism, 
thus ruling out tho possibility of a pure vacancy mechanism* 
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CHAPTER 1 


INTRODUCTION 

1.1 IMPURITY DIFFUSION AND DIFFUSED PROFILES IN SILICON 
TECHNOLOGY 

The terminal parameters of a silicon device/l.C. are 
dominantly influenced by the impurity distribution in the 
active region of the device. These impurities are introduced 
into silicon by thermal dif fusion/ion implantation, the latter 
becoming increasingly popular because of the higher control it 
permits in the impurity deposition process* 

The deposition is only one of the several high temperature ' 
processes used in the fabrication of integrated circuits, others 
being drive-in, oxidation and dielectric deposition etc. In 
the fabrication process sequence, the deposition is often 
followed by one or more high temperature processes. The doposi-J 
tod impurities redistribute themselves via thermal diffusion 
during these subsequent process steps, the extent of redistri- ; 
bution being a function of the process ambients and the tempera- 

I 

tures. The situation is further complicated by the facts that 

(i) bipolar device processing, in certain conditions, involves ; 
the simultaneous diffusion of donor and acceptor impurities, ! 

i 

which significantly interact with each other during diffusion ^ 

! 

(ii) the diffusion is a strong function of impurity concentra- 
tion, especially at high concentration values, and (iii) the 
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impurity diffusion is influenced by the lattice strain that 
may be caused by the deposited or grovm surface films or 
due to the misfit of impurity atoms in the silicon lattice 
[1,2], The impurity distribution ultimately obtained in the 
device is the net result of all the process variables mentioned 
i.e., temperature, impurity concentration, lattice strain, 
ambient etc* 

A wealth of experimental information is available in 
literature on various aspects of dopant impurity diffusion in 
silicon [l3« The basic question of the diffusion mechanism 
of dopant impurities in silicon, however, is still not answered 
uniquely. This has been and continues to be, a major limita- 
tion in developing the mathematical models for diffusion in 
silicon, which may be used to predict impurity profiles for a 
given set of process variables* 

The last seven years have witnessed a renewed interest in 
understanding and mathematically modelling the dopant impurity 
diffusion in silicon [l-7]. Efforts made in this direction are 
specifically aimed at developing comprehensive mathematical 
models, consistent with the experimental results, to accurately 
relate the process variables with the resultant impurity pro- 
file. The reasons for this renewed interest in the area are i 
(i) The trial and error method, hitherto used successfully, to 
obtain the desired impurity distribution, becomes costlier and 
more time consuming as the IC complexity increases, (ii) To 
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shorten the time required for the development of a new product, 
the computer aided design methods are to be used. This, in 
turn, necessitates the availability of accurate mathematical 
models, (iii) It is now well recognised that the reason for 
our limited capability of accurately relating the process 
variables with the dovico/circuit performance lies in in~ 
adequate process~profilo models rather than devico/circuit 
models [l]. 

The objective of this work is to study the effect of 
process ambient, i.e., presence of oxygon and chlorine, on the 
diffusion of phosphorus in silicon. Various phenomenological 
models proposed to formulate the effects of oxidation have 
also been examined. As will bo soon in Chapter IV, the 
effect of chlorine on phosphorus diffusion also yields direct 
evidence about the diffusion mechanism of phosphorus in silicon 
as suggested by Goselo and Frank [5]. The following section 
gives a brief review of the earlier work done in the area. 

1.2 PHOSPHORUS DIFFUSION IN OXYGEN AND OXY-CHLORO AMBIENTS 

1.2.1 Oxidising Ambient 

Tho offset of the prosonco of Og in diffusion ambient was 
first reported by Nicholas [9] who observed larger junction . 
depths in p-silicon samples diffused in partially oxidising 
ambient as compared to the same in samples diffused in an inorl; 
ambient. In his experiments, tho surface in^urity concentratio 
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was high ( > 10 cm ) and high concentration effects, naturally, 
were also present. Later, Ghoshtagore [lO] carried out on ela- 
borate set of experiments, diffusing phosphorus under intrinsic 
diffusion conditions and confirmed that the present of O 2 in 
diffusion ambient enhances the phosphorus diffusion coefficient 
in silicon. 

A systematic and quantitative evaluation of oxidation 
enhanced diffusion (OED) of phosphorus was reported by Masetti 
et al. [11-123* His main observations wore : (i) The phos- 
phorus OED in silicon is proportional to oxidation rate at a 
particular temperature, (ii) The proportionality constant is 
less for the caso of (O 2 + steam) oxidation than for oxidation 
in dry O 2 . (iii) The OED is temperature dependent and in- 
creases with decreasing temperature. The surface concentra- 
tion in his experiments, however, wore of the order of 
19 —3 

5x10 cm and the high concentration effects may bo expected 
to have contributed to the estimated diffusion coefficient 
values* Subsequently, Tnniguchi et al. [13] carried out more 
detailed study of OED of phosphorus (as also of boron) in 
<100> silicon. He reported that (i) at a particular tempera- 
ture OED depends only on the oxidation rate and is independent 
of oxidising agent as reported by Masetti et al., (ii) the 
OED is affected by inpurity concentration for values greater 
than approximately the intrinsic carrier concentration, (iii) 

The OED has a finite range in silicon, estimated to be about 
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25 p,m from the oxidising surface. This was an important 
observation as it indicated that the defect species i ,e, self- 
interstitials responsible for OED were able to penetrate 
only up to a certain depth. In both these experiments, the 
OED was studied by observing the redistribution of phosphorus, 
deposited by CVD, during oxidation. For estimation of diffu- 
sion coefficient, the initial distribution was approximated 
with an ideal distribution, i.e., erfc or Gaussian distribu- 
tion. Antoniadis et al. [14] and Lin et al. [15] reported 
more accurate estimates of OED of phosphorus. They studied, 
redistribution of low dose ion implanted phosphorus. For 
estimating the diffusion coefficients, the diffusion equation 
in silicon was numerically solved using the measured (rather 
than approximated) initial impurity distribution. The studies 
on orientation dependence of phosphorus OED have not been 
reported in literature, although quite a few reports on the 
orientation dependence of boron OED arc available [6]. Oxida- 
tion retarded diffusion (ORD) has also been reported by Francisi 
Dobson [17]. This effect has been observed at high tempera- | 
tures (> 1150 C) and long oxidation/dif fusion times* It may ; 
bo pointed out that for short durations (_< 1 hrs) OED is i 

observed oven at those temperatures. 1 

As indicated in Sec. 1.1 the modelling of OED would 
require the knowledge of phosphorus diffusion mechanism in 
silicon. Although the point d^f-ects responcifoip for diffusion 
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in inert ambient are still not uniquely identified, the same 
responsible for the OED have been determined to be self- 
interstitials which are generated during oxidation at the 
Si 02 -Si interface and diffuse into the bulk silicon subse- 
quently [5,6], This information is derived from the fact that 
under the same conditions in which the OED is observed, the 
oxidation induced stacking faults (OSF' s) have been found to 
grow [6], These OSF's have been determined to be extrinsic 
type which could only grow by absorption of self-interstitials. 
In an inert ambient these OSF’s have been observed to shrink. 
This indicates the presence of self-interstitials in super- 
saturation during oxidation. Various models have been proposed 
for the kinetics of these excess self-interstitials to quanti- 
tatively explain the OED as well as OSF growth data [l9-2l]. 
These will bo discussed in detail in Chapter IV in context with 
a multitude of important experimental observations on OED and 
OSF growth. However, it may be pointed out here, that kinetics 
of generation, motion and annihilation of excess self-inter- 
stitials is still not precisely known and consequently the 
modelling of OED and OSF growth is in a transition stage. 

1,2,2 Chloro-Oxidising Ambient 

The effect of chlorine on phosphorus (and also of boron) 
diffusion in oxidising ambient has been reported in detail 
only by Habqta ot al. [22] . The chlorine was observed to cause 
a reduction in OED, This observation is consistent with the 
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observed retardation in OSF growth in chloro-oxidising ambient 
[6], The reduction in OED was observed to be constant for 
HCl concentration greater than about 4/. and completely countera- 
cted the OED at temperatures 1100°C-1150°C. These data, how- 
ever, covered only a narrow temperature range, lOOO^C-llSO^C 

20 —3 

and involved high impurity concentrations, i.e., 2 • 

Their observations included an anamalous feature, that is, the 
larger diffusion coefficient of phosphorus in <111> than in 
<100> Si in oxidising ambient at lOOO^C. This has been attri- 
buted to the extrinsic diffusion conditions of the experiment 
[18]. Further experiments to evaluate the effect of chlorine 
on phosphorus diffusion, free from other effects, are obviously 
called for. 

1.3 THIS WORK 

As mentioned in Sec. 1.2, the experimental results on OED 
reported in literature, have been obtained by studying the 
redistribution of phosphorus during oxidation. The extraction 
of diffusion coefficient from the redistributed profiles is 
rather an involved process because this experimental situation 
is not amenable to analytical mathematical methods. Further- 
more, the approximations made for the initial distribution or 
errors in the experimentally determined initial distribution are 
likely to introduce errors in the estimation of diffusion 
coefficient. 
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An alternative method to study the effect of ambient 
on phosphorus diffusion is to carry out diffusion from a gaseous 
source, in oxidising/chloro-oxidising ambient. Mathematical 
analysis of this situation is very sinplo. Since phosphorus 
tends to segregate in silicon, it is possible to carry out con- 
stant source diffusion in very highly oxidising ambient 
( 9 ^°/. O2), thus permitting the observation of the OED of the 
same magnitude as may be observed under dry O2 oxidation condi- 
tions. The simplicity of the theoretical treatment of this 
experimental situation permits an easy determination of diffu- 
sion parameters. 

This work utilizes the above described technique to study 
the phosphorus diffusion in oxidising and chloro-oxidising 
ambient. The plan of the work is as follows : 

In Chapter II, a phenomenalogical model for the diffusion 
in oxidising ambient is presented. The associated moving boun- 
dary problem has been solved for impurity distribution in sili- 
con, At the PSG/Si interface the impurity concentrations in 
PSG and Si are assumed to approach thermodynamic equilibrium 
gradually rather than instantaneously. 

The details of the experimental work are presented in 
Chapter III. A description of the experimental set up, plan of 
diffusion and measurements carried out on the diffused samples, 
is given. The procedures adopted for the determination of 
impurity distribution in silicon and extraction of diffusion 
parameters are described. 
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The experimental results are discussed in Chapter IV. 

The special features of the diffusion in oxidising ambient are 
pointed out. Various models proposed for OED and OSF growth 
have been examined in the light of the observations made in this 
work. The observed effects of chlorine on phosphorus diffusion 
are discussed and it is pointed out that the lowering of diffu- 
sion coefficient in chloro-^oxidising ambient below its value 
in inert ambient negates a pure vacancy diffusion mechanism of 
phosphorus in silicon. 

Lastly, Chapter V presents a brief summary of the work 
presented in the foregoing chapters. 
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CHAPTER II 

DIFFUSION IN OXIDISING A^ffiIENT ~ 

THEORETICAL MODEL 

Mathematical formulation of the process of impurity 
diffusion from a constant source, in oxidising ambient, is 
exactly similar to the problem of impurity redistribution 
in a homogeneously doped sample during thermal oxidation. 

The solution to the latter has been proposed by Grove et al. 

[l]. In their model, however, the impurity concentrations in 
Si02 and Si at the interface, have been assumed to be in che- 
mical equilibrium, right from t = o . This assumption, natura- 
lly, leads to the prediction of a time-invariant impurity con- 
centration at silicon surface. In case of boron redistribution 
during steam oxidation, the surface concentration is reported 
to be time invariant, in agreement with the prediction of 
Grove’s model [2]. For the case of phosphorus [3] and boron 
[ 4 ] diffusion in weak oxidising ambient, however, the surface 
concentration has been reported to be time-variant, approaching 
a constant value for long diffusion times. This indicates that 
under slow oxidation conditions, the impurity concentration at 
silicon surface approaches its equilibrium value gradually 
rather than instantaneously- No experimental results appear to 
have been reported in between these two extreme cases. It is 
however evident that in general the inpuritios in SiOg and Si 
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at the interface approach thermodynamic equilibrium gradually, 
the rate of approach to the equilibrium being a strong function 
of oxidation rate. 

In the Section 2.1, the mathematical formulation of impu- 
rity diffusion in oxidising ambient is presented. A rate 
limited flux boundary condition has been assumed at the inter- 
face, to model the gradual approach to equilibrium, as first 
proposed by Smits and Miller [S] . In Sec. 2.2, the solution to 
problem posed in Sec. 2.1 is discussed. 

2.1 MATHEMATICAL FORMULATION 

The process of impurity diffusion during thermal oxidation 
is illustrated in Fig. 2.1, Using a two coordinate system, 
shown in the figure, the diffusion eqs. in dopant-silicate 
glass (henceforward referred to as Si 02 in this chapter) and 
Si can be written as, 


dc,{x,t) d^C,(x,,t) 

! Xix„(t) 


( 2 . 1 ) 


and 


dC2(y,t) fii^C^Cy,!) 

= D — _£ , r,X^(t) < y < oo (2.2) 

' t > o 




Here oxide growth has been assumed to be parabolic so that 
X^(t) == The boundary and initial conditions for the above 

equations for the case of diffusion from a constant dopant 
source are given below. 



Impurity Cone- 


12 



~Cb 


FIG* 24 Phosphorus Diffusion in Silicon in 
Oxidising Ambient . 
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(i) At Si02 surface, C 2 _(o,t) - = const. (2,3) 

(ii) Rate limitation at Si-Si02 interface : 

Rate limitation at the interface is expressed as, 
dC2 

[-□2 r .X^Ct) ^I^’^‘^l»x=XQ(t)‘"^2ly= r .X^Ct)^ 

(2.4) 

or 

dC^ ■ 

“^2^y= x^XQ(t) “ (2.4a) 

whe re , 

a = K/D 2 
and f(t) = 

(iii) The condition of impurity conservation in Si-Si02 is 
given by [ l], 

dCi dX dC2 dX 

'^l]x=Xg(t) = ■57- •'^^y=r.Xo(t) 



( 2 . 5 ) 

silicon, 


C2(”,t) = -Cg 

( 2 . 6 ) 

{ y f 0 ) “ 

( 2 . 7 ) 


2.2 IMPURITY DISTRIBUTION IN SILICON 

A general solution to eqn, (2,1), for parabolic oxide 


growth, may be written as [6], 
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C,(x,t) = S (Dj^t)"/^ [a h (— ^ 
n=o ^ " " 2fD^t 


) + erfc (- 


■)] 


( 2 , 8 ) 

The coefficients may be eliminated by using eqn. (2.3). 
The resulting expression for Cj^(x,t) is, [Appendix l], 

C,(x,t) = a +(C -a ) erfc( — ~ — ) + Z (D^t)*^^^ a h (-^^ — ) 

X U U U rN*vrr^ x onfrv X 


2V'Dj_t n=l 

Z [ni .2’^] i’^ erf c(-~3^^-~) 


2fDTt 


n=2 


2fDjLt 

(2,9) 

From eqn. (2.9), f(t) of eqn. (2,4) and L.H.S. of eqn. 
(2,5) may be written as, 


f ( t) = -maCj_(xQ,t) 


-raa[ aQ+(CQ— S q) erf c( 6 ) + Z (Dj^t) ^ a^^ hj^( 5) 

n=l 

(Dj^t)’^/^- (2’^.nOa^.i’^ erfc( 6 ) 


n=2 

where 6 == Y’B/(2V’Dj^) 
and, 

dC, dX 


(2.10) 

( 2 . 11 ) 


fD 


[Dl ^ ^ aQ+(CQ-ao){26 erfc( 6 V 

+ Z (Dit)"^/^ ®nVl^*^ (Dj^t)’^/^ 2^.nl 

n=l n=2 


2e 




2(n+l)i'^'^^ erfc(6) 


( 2 . 12 ) 



15 


To determine C 2 (y,t), the impurity distribution in 
silicon, eqn. (2.2) is first transformed to (z,t) domain 
where z = y'-rX^Ct). The transformed diffusion eqn. and the 
related boundary and initial conditions, after substituting 
for X^Ct), may be written as, 


^2 


d^C2 rfB 
dz^ 2ft 




(2.2a) 


dC^ 

= -m.a.Ci(XQ(t),t) = f(t) (2.4b) 


C2(<«,t) = -Cg (2.6) 

^2^^^^^ ~ (2.7a) 


The general solution to eqn. (2.2a) with boundary and 
initial conditions (2.4b), (2»6a) and (2.7a) is [see 
Appendix II] 


€ 2 ( 2 , t) - -Cg + f g( z, t-X) .p( X) .dX 

o 


(2.13) 


where, 

g(z,t) 


„ r . .z+rfBt^ a(z+rfBt)+a%t z+rfBt 

p . [ e r f c ( ) -«e .erfc( 


+ afD2t)] 


2fD2t 

(2.13a) 


P 




(2.13b) 


er 


p(t) = +Cg + i f(t) « +m C 3 ^(X^(t),t) 


(2.13c) 
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Using eqn. (2.13), the expression for R.H.S. of eqn. (2.5) 
is determined (Appendix II). 

^ + IB ^2^=0 = m + *=8 “Cv 

+ (CQ-a^) erfc(6)]} . { (2D2ft+ -^1^) .S( t)+rirB} 

+ E m h^(6) £(2D2ft+ S(t~X)dX 

n=l o 

« 

+ rVB.W^ } - ^e"° m a 2“^ nl i’^ erfc(6) £ ( 2 D 2 Vt 

n=2 

+ ■^) It S(t-x)dX + rfBt’^/^ } ] (2.14) 

Now eqns. (2.12) and (2.14) may be combined in eqn. (2.5) 
to write the impurity conservation condition in terms of the 
coefficients which are the only unknown quantities in 

eqn. (2,13), the equation for impurity distribution in silicon 
(see Appendix III for eqn. (2.5)). Because of the complicated 
nature of the equation, it is not possible to determine Sj^'s 
for a general case from the conservation condition. However, 
for some special cases it is possible to evaluate a^’s. These 
are considered below. 

Case I ; Negligible PSG Growth, 6 is o 

For this case, the expression for p(t) is determined as, 
(see Appendix III). 

p( t) - m Cq 


(2.15) 
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After substituting B = o in eqn. (2.13), the expression for 


impurity distribution in silicon can be written as. 


C2(z»t) = -Cg+mCQ[erfc( 


2VD2t 


az+a'^D^t 

■)-e erfc( 


2fD2t 


aVD2t)] 


(2.16) 


Equation (2.16) is identical to the one reported by Smits 
[S] for diffusion from a constant vapour source in inert ambient. 
It is an expected result because as 6 *• o, the experimental 
situation approaches that of diffusion in inert ambient from a 
constant source. 


Case II : Large Value of Rate Constant/Diffusion Time,a'fD2t >> 1 

The expressions for p(t) and g(t) are derived to be, (see 
Appendix III) 


p(t) = A = m[aQ + (C^-a^) erfc(6)] 

,2 — r^B/4Do 

m[CQ(e~ ,e .erf (6).Cg] 

= ^ 

( e~ /Y‘7t)+ 6 .erf (6)+m.r.6 .erfi^(6) .F(r,B,D 2 ) 

(2.17) 


whore, F{r,B,D2) = -1 + 


1 


f7t(rfB/2fD2).e^ .eTfc{xtB/2fD^) 


and 


g(t) « p[erfc( 


rlTB 


•) 


2fD^t 2fQ^ afuD^t 

mum mm m mmm 


] 


(2.19) 
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Substituting for p(t) and g(t) in eqn. (2.13), the 
expression for C 2 (zjt) is written as, 

C^(z,t) = -Cr,+A[erfc( — ~ + — I ] (2.20) 

“ 2VD2t 2fD2 (xfnD2t 

In eqn, (2,20), if a -* «, it is seen that, 

C^(z,t) = -Co+A erfc( — - — + ^^) (2.21) 

^ 2fD2t 2^02*'' 

Equation (2,21) is exactly same as reported by Grove et al. 
[l] for redistribution of impurities in homogeneously doped 
sample. This is because as a •* oo, Cj^ (C 2 /m) at the inter- 
face in order to maintain a finite impurity flux across the 
interface. The relation Cj_ = C 2 /m essentially represents the 
condition of thermodynamic equilibrium at the interface, assumed 
in Grove's model. 

Case III : otfDot » ^ 

2fD2 

Since at the diffusion temperatures normally encountered, 
rfB « aiTDj, this condition represents a practical diffusion 
process. The expression for p(t) is derived as, (see 
Appendix III) 

h„(6) „ « «/o 

p(t) sa A + E -cp — r \ 2m Apf! n (-1)^ ( 0 :^ 02 )" t^'^ 
n=l n+l'‘°'' 


( 2 . 22 ) 
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where 


2 

C^Ce ^ 6 .erf(6) Cg 

( e“ ^ /V^t) + 6 .erf ( 6 ) 


(2,22a) 


Using equations (2.22) and (2.13a), the expression for 
C 2 (z>t) can be written as, 


C 2 ( z,t) 


-Cg+A.g( z,t) 


+ I J[1 


2mA\xp n (-1)’'^ . 


(oTDj)"] |t / g( 2 ,t-X )dx (2.23) 

It can be seen that series in eqn. (2.23) converges only 

for small values of aY’D 2 t such that aY’D 2 t < 1. Since the value 

3 4-1 

of 'a' in weak oxidising ambient is typically 6x10 -10 cm 
[7], this equation is useful only for small diffusion lengths 
(less than a micron). 


The three cases considered above, do not encompass all 
the possible diffusion situations. Case III, however, corres- 
ponds to an important situation namely deposition process which 
generally involves small diffusion lengths. In the following 
chapters, we will be concerned with Cases I and II only. 
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CHAPTER III 

PHOSPHORUS DIFFUSION AND CHARACTERIZATION OF DIFFUSED 

LAYERS 

In the theoretical model presented in the preceding 
chapter, a constant dopant source and position independent 
diffusion coefficient were assumed. These assumptions impose 
certain constraints on the experiment, which have been discussed 
in Sec. 3.1. In Section 3.2, the details of experimental set- 
up and the scheme of experimental work have been given. 

Sections 3.3 and 3.4 are devoted to the characterization of PSG 
and the phosphorus diffused layer in silicon, respectively. The 
last section i.e. Sec. 3.5 describes the method of determinatiorl 
of diffusion parameters. 

3.1 EXPERIMENTAL CONSTRAINTS 

In order to closely match the experimental situation with 
the theoretical model proposed in Chapter II, the following | 
conditions should be satisfied. 

(i) The partial pressure of the dopant in diffusion ambient | 
should remain constant with time. This requires that the 
dopant should have good miscibility with the carrier/dilutent 
gases (N2»Ar,02 etc*). Furthermore, the deposition of dopant ; 
impurity on the walls of diffusion tube at the cool input end 
should be minimised as it would act as a time dependent second-! 
ary dopant source. 
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(ii) The dopant concentration at silicon surface should be 
less than the intrinsic carrier concentration at diffusion tem- 
perature. It ensures that the diffusion coefficient is impurity 
concentration (and therefore position) independent as assumed 
in the model [ l] . 

A reasonable choice of dopant source satisfying all the 
above requirements is P [3]. The experimental set-up used 
for diffusing phosphorus using this source is described jn the 
following section. 

3.2 PHOSPHORUS DIFFUSION 

3.2.1 Experimental Set-up 

A schematic of the experimental set-up is shown in Fig. 3.1* 
The diffusion furnace used is Thermco Mini-brute type MB-71. As 
per manufacturer's claim, this furnace provides a flat tempera- 
ture zone of 12-14” with spatial variation within + 0.5°C, with 
a quartz tube of OD 3,25”. The temperature in the flat zone is 
expected to be stable within + 0.25°C. In the experiments 
reported in this work a quartz tube of 1,5" ID was used and 
wafer carrier was a quartz boat of 3" length. The diffusion 
temperature was measured using a Pt.~Pt. + 10:4 Rhodium thermo- 
couple. The quartz tube was so positioned in the furnace that 
the gas entrance point (see Fig. 3.1) is maintained at a teirpe- 
rature > 500°C. This eliminated the deposition of P 2 O 5 at the 
input end to the extent that even after 15 hours of diffusion 



/Furnace 
r— Liner 



FI6.34 Schematic of experimental set-up 

for Phosphorus diffusion using PH 3 
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there was no visible layer of P2’^5* However, the ^ 2^5 deposi- 
tion was observed at cooler output end. 

Teflon tubes were used for making gas connections. The 
dopant source was an electronic grade PH^+Ar mixture supplied 
by Matheson/international. Other gases used in diffusion 

were also of high purity, having a maximum impurity 
concentration of less than 6 ppm. The needle valves used to 
control the flow of gases were made of teflon. 

The experimental set-up shown in Fig. 3.1 corresponds to 
the case of diffusion in inert (i.e., negligible yi of O 2 ) 
ambient. For diffusion in strong oxidising ambient, the H2»02 
gases lines were interchanged at the input end of the res- 
pective flowmeters, N 2 flow was blocked and mixer flask was 
bypassed. For diffusion in chloro-oxidising ambient, the 
trichloroethylene (electronic grade) was filled in the mixer 
flask and N 2 was bubbled through it. One inlet of mixer flask 
(bubbler) was blocked. At the outlet of the bubbler a two 
way teflon valve was connected so that the TCE carrying N 2 
could either be exited to atmosphere or let into the input gas 
line of the diffusion tube. 

3.2,2 Details of Experiments I 

The dopant source used in the experiments was electronic 
grade PH^+Ar mixture (1000 ppm PH^ in Ar) . The diffusions were' 
carried out in three different ambients. 
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(i) Strong oxidising ambient : 90.1 0^ + 9.9 Vyi (PHg+Ar) 

(ii) Inert ambient : 95.5 N 2 + 1.5 O 2 + 3.0 V><(PH 0 +Ar) 

(iii) Chloro-oxidising ambient ; 94.5 O 2 +O .5 V>< (PH^+Ar) + 

5.0 N 2 bubbled through trichloro- 
ethylene at 27°C. 

The ambient compositions given above are the typical ones. 
The exact ambient composition for each sample is given in 
Tables 3.1 - 3,4. The pre-diffusion cleaning of silicon sam- 
ples was done using hydrogen peroxide based cleaning solutions 
[ 4 ]. Prior to loading the samples into the furnace, the gases 
were allowed to flow through the tube at the p re-determined rate 
to let the flow of respective gases stabilize. The quartz boat, 

with samples horizontally placed on it, was quickly loaded into 

* 

the furnace without interrupting the gases flow. However, the 
dopant gas flow was stopped at the time of unloading. 

The samples used in these investigations were squares of 
12 mm X 12 mm size. Starting wafers for <111> samples were p- 
silicon slices of 3" dia. and 5-10 Q .cm resistivity. The same 
for <100> samples were p-silicon wafers of 2*’ dia. and 4-6 Q.cm/| 
22-50 Q. cm resistivity. The starting wafers were supplied with 
one side mirror polished, ' 

„ 

In case of diffusion in the inert ambient, the dopant gas was 
interrupted while loading the sample. 
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A typical sample load for a diffusion consisted of two 
<111> and two <1CX)> samples of the type mentioned above, along 
with a small size (typical dimensions 5 mm x 3 ram) p-sample, 
ear marked for oxide thickness measurement. In case of diffu- 
sions in inert ambient, this last sample was omitted. 

Tables 3.1 - 3.4 give the details of each diffusion, such 
as diffusion temperature, time and important sample specifica- 
tions along with sample identification numbers. First the 
diffusions were carried out in strongly oxidising ambient. 
Appropriate diffusion times were chosen to allow the growth of 
2000°A - 3000°A PSG (see Table 3.1). The rate of PSG growth 
in these diffusions was found to be only slightly lower than 
the oxidation rate of silicon in dry O2 at corresponding tem- 
peratures. The second set of diffusions were done in inert 
ambient for the same time interval as chosen in case of strongly 
oxidising ambient at respective temperatures. The presence of 
the PSG film in these cases was visually undetectable. Since 
the oxidation rate of silicon is proportional to the partial 
pressure of oxygen in the ambient, tho PSG film thickness in 
these cases is estimated to be less than lOO^A. The difference 
in phosphorus diffusion coefficient in above two cases, there- 
fore, can be attributed to the oxidation of silicon during 
diffusion in strong oxidising ambient. 

The next two sets of experiments wore performed to investi- 
gate the effect of chlorine on phosphorus diffusion (see Tables 
3.3 - 3.4). The two sets of diffusions wore done in chloro- 



Table 3.1 Diffusions in Strong Oxidising Ambient Long Diffusion Times 
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used was ordinary grade and is expected to contain 5^ N2 


Table 3.4 Diffusions in chloro-oxidising Ambient 
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oxidising and strong oxidising ambient respectively. The 
diffusion times for these experiments were selected to be 
equal at respective temperatures and considerably shorter 
than the same for the diffusions described in preceding 
paragraph. The effect of chlorine on phosphorus diffusion 
in silicon can be obtained by comparing the diffusion coeffi- 
cients estimated in the two cases. Since for strong oxidising 
ambient case, the diffusions have been done for two different 
time intervals, a comparison of respective diffusion coeffi- 
cients also allows the determination of time (and hence oxida- 
tion rate) dependence of OED. The shorter diffusion times of 
these experiments result in reduced junction depth which in 
turn also reduces the time taken in evaluating the impurity 
distribution in diffused layers. 

3.3 CHARACTERIZATION OF PSG 

The PSG parameters of interest are C^ and X^{t) . The 
knowledge of X^Ct) is essential for determination of D2. As 
will be explained in Chapter 4 (see Sec. 4.1), the knowledge 
of Cq provides a way to calculate the value of D^. The experi 
mental determination of C^ and XQ(t) is described below. 

3.3.1 Estimation of C^ 

0 

Approximate value of Cq was determined from the etch 
rate of PSG in p-etch [5,6]. The calibration curve given in 
ref- was used to convert etch rate of PSG in p-etch to 
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mole percent of P 2‘^5 PSG. The average value of etch rate 
of PSG was found to be 4.3 + 0.2 A/sec. for all the cases of 
diffusions in strong oxidising/chloro-oxidising ambient. The 
corresponding concentration of P 2 O 5 PSG is 4.0 mole percent. 
If the density of PSG is taken to be same as that of Si 02 r 
thermally grown in dry O 2 (i.e., 2.25 gm/c.c), the concentra- 
tion of P 2 O 5 iri PSG - and the value of Cq are calculated to be 
1.7x10*' cm and 3.4x10*^ cm respectively. 

It would appear that the value of determined as above, 
is the average value of phosphorus concentration in PSG rather 
than the phosphorus concentration at the PSG surface- It has 
been shown, however, by radiochemical analysis, that the phos- 
phorus distribution in PSG is box type [7] and for a box type 
distribution, the average value of phosphorus concentration will 
be nearly equal to that at the surface. The above determined 
value of Cq is therefore, expected to bo a good approximation to 

true value of C^. 

o 

3.3.2 Determination of PSG Thickness 

Ellipsometry, MOS capacitance measurements and interfero- 
metry are the most commonly used methods for measuring the oxide 
thickness. In this work, the oxide thickness was measured using 
interferometry. 

The basic principle of interferometry, as applied to the 
measurement of step height is illustrated in Fig, 3.2. The 
specimen is covered with a reflective coating of silver and a 



Monochromatic tight 


FIG. 



(b) 


.2 Oxide Thickness Measurement by 
Interferometry (a) A typical position 
of reference plane with respect to 
specimen ( b) Inter fcrcnce fringes 
as seen through microscope. 



33 


partially reflecting plane is positioned over the specimen 
as shown in Fig. 3.2a. V/hen a monochromatic light falls on 
the assembly, an interference fringe pattern, as shown in 
Fig. 3.2b, is created due to the path difference of the lights 
reflected from the reference plane and the specimen surface. 
These fringes are viewed through a microscope. The step height 
tg is related to the fringe displacement D, the normal fringe 
spacing d and the wavelength of the monochromatic light Xo, by 
the relation, 



The measurements were carried out on'^A scope model 960-4020 
(Varian make). 

3.4 Characterization of Diffused layers 

Differential sheet resistance method has been used to 
evaluate the diffused impurity profiles. The method is appli- 
cable only to the cases in which majority carrier concentration 
distribution is same as impurity distribution. In case of phos 

phorus in silicon, the two distributions coincide for impurity 

20 —3 

concentrations _< 10 cm [IJ. The impurity concentrations 
encountered in this work, are more than one order of magnitude 
lower than this limit, 

3.4.1 The Differential Sheet Resistivity Method-Theory 

The bulk conductivity at a depth ’x’ from the surface is 
related to the sheet conductivity at the same depth by the 
equation [9] 
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X . 

a (x) = S cy(x) dx (3.2) 

® x 

Differentiating eqn. (3.2) w.r.t. ’x', 



or 

G.^i {n(x)} . n(x) = - [in Pg(x)] (3.3) 

Due to complicated nature of tho function p[n(x)], eqn. 
(3.3) can not be solved analytically for n(x) and therefore 
has to be solved numerically [Sj. 

It is evident from eqn. (3.3) that tho sheet resistivity 
is required to bo measured at different depths from the surface 
by successive removal of thin layers of silicon. Tho accuracy 
of the method, therefore, depends upon (i) accuracy of the 
measurement of sheet resistivity and (ii) uniform removal and 
accurate determination of, the thickness of silicon layer * 
removed before each measurement. 

Two methods, namely, four point probe and Vander Pauw 
methods have been used for measurement of sheet resistivity in 
this work. 

(a) Four Point Probe Method 

Shoot resistivity in this case is given by [s] , 

p V c(D /s).F(x./.s) (3.4) 

® ln2 ^ ® ^ 

where I is the current flowing through tho outer probes and 'V’ 
the corresponding potential difference created between the two 



35 


inner probes by the current I,C(Dg/s) and F(Xj/s) are corre- 
ction factors. C(D /s ) takes into account the finite lateral 
dimension of the sample and is a function of sample diameter 
'D ’ and probe spacing 's'. F(x_-/s) is the correction factor 

s J 

to be included v\tien sample thickness x^ is large compared to 
's' . For Xj/s < 0.4, F is unity. The values of C and F as a 
function ofDg/s and x^/s are tabulated in ref. [9] and the 
same have been used in this work* 

To measure the sheet resistivity, a four point probe 
having a probe spacing of 25 mils (supplied by Kulicke and 
Soffa) was used. For the sample shown in Fig. 3.4a, the eqn. 
(3.4) reduces to, 

Pg = 4.1716 ^ (3.5) 


(b) Vandor Pauw Method ; 

Referring to fig. (3.4b), tho shoot resistivity is given by 


[ 10 ], 


‘’s f%2,34 f^3,14> ^<15^ 


12.34 


1 14 


) 


(3.6) 


R 


= !34 
12,34 ^12 


R, 


V, 

23,14 “ T 


(3.6a) 


23 

14 


where is potential difference between terminals 3,4 when a 
constant current lj^2 flows between terminals 1,2 and V23 is the 
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p.d. between terminal s 2,3 for current flowing through 
terminals 1,4. The c*orrection factor f is given by the trans- 
cedental equation, 


l4-exp[ ^ 


( ^12,-34 ^23.14 
% 2 ,: 34 '^^ 23,14 


)] = exp[ 


2 In 2 



R 

^12,34 ■ 

% 2 , 34 '‘‘^ 3,14 


(3.7) 


Equation (3.6) iss valid if the point contacts are placed 
at the periphery of tihe sample to make current and voltage 
measurements. All priactical contacts, however, have finite 
dimensions. In Fig -:3.4b, considering the heavily doped 
region as the contact region, the sum of the contact lengths 
of all the four con-ta»cts of the sample, is less than 5^ of 
the sample poripher'^. For these samples, therefore, no correc- 
tion needs be applied for finite contact size [lO]. 

3.4.2 Detemination nof Sheet Resistivity Distribution 


As mentionod inSSoc. 3.2.2, the phosphorus diffused p- 
samples are squares orf size 12 mm x 12 mm. For the measure- 
ments, the mesas of c*onf igurations shown in Fig. 3.4 were 
etched on those samplaes. The sheet resistivity measurements, 
each followed by renrion/al of thin silicon layers by anodic oxi- 
dation, were dome on -these mesas. 

(a) Sample Prepara-tioon : 

The process sequeence leading to the samples shown in 
Fig. 3.4 is given in -the flow chart of Fig. 3.3. The flow chart 




High Cone. Sftlecii! Vi 
Phosphorus 
Diffusion for 
Making ohmic 
Contocts 


Metal Deposition 
for ohmic contacts 



FIG. 3.3 Process Flow Chart for 
Preparation of Vander 
pauw and Four Pt. Probe 
Samples* 

















FJG.3.4 (a) Four Pt. Probe Sample 

(1)Ele vation 
(■|i)Top view 

(b) Vander Pauw Sample -Top view 




describes the process sequence used to fabricete Vander Pauw 
Samples (VDP). The processes used for preparing four point 
probe samples (FPP) are merely a subset of VDP process soquonce 
and have been identified by double-lining the blocks, each 
representing an individual process, in the flow chart. The 
essential details of the various processes, labelled as a,b,c,.. 
etc. in the flow chart, are given below. 

(a) Sample cleaning was done using the hydrogen peroxide based 
solvents described in reference (4). 

(b) The sample was then thermally oxidised in dry O 2 at 1200°C 
for 1 hr., to grow about 2000‘^A thick Si02. This oxide acts as 
a mask against the phosphorus diffusion during the selective 
high concentration phosphorus diffusion that was carried out in 
stop (d) to form ohmic contact regions. 

(c) Windows were etched in SiO^ grown in (b), using photo- 
lithography, to expose the silicon surface where high concentra- 
tion phosphorus diffusion was to bo performed, for making ohmic 
contacts. 

(d) High concentration phosphorus diffusion for ohmic contacts 
was performed using phosphosilica film (10 cm ) as a source.' 
The source was applied only to exposed silicon surface of ohmic 
contact region. The diffusion was done at 875°C for 15 min. 

I 

(o) The masking oxide and residue of phosphosilica film wore 
removed and sample cleaned by process (a). 
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(f) Now the phosphorus diffusion using PH^ was carried out 
as per details given in Tables 3.1 - 3.4. 

(g) The windows were etched in the PSG, grown in Step (f), to 
deposit metal contacts, the windows were quarter circles of 
1.5 mm radius, symmetrically placed at the four corners of the 
sample as shown in Fig. 3.4b. The PSG on the back side of 
sample was preserved during window etching. 

( h) The electroless nickel plating was done to deposit metal 
for making ohmic contacts [ll]. Since the whole sample surface, 
except the window regions, was covered with PSG, nickel was 
plated only over the bare silicon surface of window region. 

(i) This step consisted of Cr-Au evaporation at a pressure of 
exlO""^ torr with sanple heated to 275°C. This process had to 
be added because Cr-Au layer provides extra protection to the 
underlying PSG and Si surface and ensures the elimination of 
pin holes during mesa etching. 

In case of VDP samples the nickel contacts wore masked 
during Cr-Au evaporation. Heating of the sample to 275°C 
improves the adhesion of nickel contact. Furthermore, since the 
nickel was heated in the vacuum, it did not oxidise and wires 
could be easily soldered to these contacts. It may bo mentioned 
that the evaporated Cr-Au was not found suitable for making 
ohmic contacts presumably duo to higher than desirable pressure 
during evaporation. 
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For the samples diffused in inert ambient, this step was 
omitted. All such samples were evaluated by four pt. probe, 

(j) The photoresist pattern of the shapes shown in Fig. 3.4 
were generated by photolithography. The Cr-Au from the 
exposed surface was removed using etchants given in Appendix IV. 
PSG from the exposed surface was also removed by buffered HF. 

(k) Now the mesa is etched using Wright's etch [12]. This 
etchant was preferred because it gives a better finish of 
etched surface and because it has low etch rate which allows 
the etching to be done at room temperature without damaging 
the photoresist. The sample was quenched in buffered HF 
immediately after removal from etch bath and rinsed in DI water. 
The etched surface was then tested by hot probe for conducti- 
vity type . 

(l) The photoresist and Cr-Au were removed from the mesa 
region. The phosphorus diffused musa region at this stage is 
protected by PSG grown in step (f). 

(m) In case of VDP samples, thin copper wires wore soldered 
to the nickel contacts using Pb-Sn solder. 

Before measurements, PSG was removed from mesa region. 

(B) Sheet Resistivity Measurements 

The sheet resistivity measurements were carried out using 
a Kithley digital multimeter (model 191) as a Voltmeter and a 
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Kithley digital electrometer (model 616) as an ammeter. All 
the measurements were taken in the current range 1-100 ii.A. The 
measured voltages corresponding to these currents were less 
than 15 mV, thus eliminating the possibility of error caused 
due to forward biasing of underlying p-n junction [l3]. 

All the circuit connections were made using sheilded wires. 
The sample was covered with a metal box during the measurement, 
to ensure that the measurements are free from any error caused 
by photogenerated carriers. 

(C) Removal of Thin Layers of Silicon 

Thin silicon layers were removed, preceding each sheet 
resistivity measurements by anodic oxidation [14-17]. The 
anodic oxidation cell used for the purpose is shown in 
Fig. 3.5a. The eloctrolytu used was 0.08N solution of 
in ethylene glycol + 4>< water. The advantage of keeping the 
silicon sample in face up position is that the gas bubbles 
released during the oxidation process move up towards the 
surface of electrolyte and silicon surface is always in con- 
tact with the electrolyte. This leads to uniform and repro- 
ducible oxide growth. For providing the back contact to 
silicon, mercury has been used as it was found to provide 
good electrical contact to a lapped p-silicon surface of 
resistivity upto 50 cm. Because of the hydrophobic nature 
of teflon and silicon surface, the electrolyte does not creep 
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out and anodic oxidation is strictly confined to the area 
bounded by teflon structure containing the electrolyte. 

In order to ensure that silicon layers of equal thickness 
are removed in each oxidation step, a standard anodic oxidatior 
process was defined in the following manner. ’The cell voltage 
is first increased so that the anode current rises to 8 mA. 

As oxidation proceeds, the current decays. The cell voltage 
is now gradually increased keeping the current between 6 mA - 
8 mA, till the applied voltage reaches 120 V. The current is 
allowed to decay to 1.2 mA. The oxidation is then terminated 
by bringing the applied voltage to zero' . This standard pro- 
cess was used for removal of thin silicon layers throughout 
this work. 

To accurately determine the thickness of silicon removed 
in each standard anodic oxidation process (SAOP) the following 
procedure was adopted. A phosphorus diffused p-Si (12 mm x 
12 mm) sample was thermally oxidised to grow an oxide layer of 
about 26CX)°A. The Si02 was then removed from half of sample 
by photolithography. Now a SAOP was performed. The anodic 
oxidation occurred only over the exposed Si-surface. This 
anodic oxide was etched in p-etch [5] for 20 sec. Since the 
etch rate of anodic oxide is about 50°A/sec as compared to 
2'^A/sec of thermal oxide, an etching of 20 sec in p-etch was 
enough to ensure complete removal of anodic oxide while reduc- 
ing the thickness of masking thermal oxide only by about 40°A, 
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The above process of anodic oxidation was repeated 10 times. 
Even at the end of 10th anodic oxidation there was enough 
thermal oxide to mask the underlying silicon surface from 
anodic oxidation. After last anodic oxidation, the sample 
surface appeared, as shown in Fig. 3.5b(iv), A rectangular 
groove of 2 mm x 5 mm was then etched in the oxide such that 
the boundary of anodic and thermal oxides falls in the centre 
of the groove. This resulted in the exposure of silicon stop, 
height of which is equal to the silicon thickness removed in 
10 SAOP’s and also a step the height of which is equal to the 
thickness of anodic oxide grown in one SAOP. Now a silver 
layer of about 2000^A thickness was evaporated on the sample 
and the step heights measured with Vartan X Scope model 
960-4020. 

The above experiment was carried out on two test samples. 
The average silicon step height was determined to be (2440t5/X 
and average anodic oxide thickness to be 543‘^A+15°A. The sili' 
con removed per SAOP was thus taken to be 244°A and the ratio 
of the volume of silicon consumed to the volume of anodic 
oxide grown is determined to be 0.449, which is in agreement 
with the reported values of 0.38-0.45 [8,16], within experi- 
mental errors. 

3.4*3 Determination of Impurity Distribution 

Impurity distribution from sheet resistivity distribution 
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was determined using eqn. (3.3). For M-(n), the following 
expression was used [18], 


l^(n) 


55.24 + 


1388.57-55,24 


1 + (- 


n 


1.072x10 


17 


0.733 


(3.8) 


The measured P_(x) vs x data were first smoothed using 

s 

cubic spline approximation. An IMSL subroutine (International 

Mathematical Subroutine Library) ’ICSSCU' was used for the 

purpose. The extent of data smoothing was adjusted to obtain 

2 

a reasonably smooth impurity distribution. The value of R 

^ 2 

for smoothed data was between 99.9 - 99.5 . The value of R 
was 99.9 for inert ambient data, 99.8 for strong oxidising 
ambiont(0/0A - samples) and 99.5 for other cases. 

The slope d P /dx was determined using another IMSL sub- 
routine ’DCSEVU*. The impurity distribution was then deter- 
mined by solving eqn. (3.3) for n(x) numerically by Newton- 
Rap hson method. 


* 2 

The R is defined as 


2 

R = 1 


n 

E e 
m=:l 


m 


n _ , 

2 [Y -Y]‘ 
■* 

m=x 


where 


Yjjj = Observed value of dependent parameter 

y = Mean value of 

e = Y - Y 
m sm m 

Y„„ = Smoothed value of Y„ 
sm m 
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The impurity distribution and the measured sheet resisti- 
vity distributions for different samples are plotted in Figs. 

3-7 - 3.24a. 

3.5 DETERMINATION OF DIFFUSION PARAMETERS 

The values of diffusion parameters, i.e., A,D 2 and a, were 
determined in the following manner. 

First the value of a was estimated from the surface concen- 
tration vs. diffusion time data plotted in Fig. 3.6. Since it 
has been shown that the ’a’ is temperature and crystal orienta- 
tion independent [18-19], the value of a was determined only 
at one diffusion temperature for each diffusion ambient. 

It can be seen that in case of diffusion in strong oxidi- 
sing ambient, the surface concentration is time-invariant. The 
value of a for these cases, therefore, is taken to be infinity. 
In inert ambient, the surface concen-tration varies with time. 

The value of a in this case is determined by using oqn. 2.16 
at X = o. The diffusion coefficient of phosphorus in silicon 
at IISO^C was taken from Ghoshtagore* s data [19], for deter- 
mination of a. The calculated value of 'a* for this case is 
3 —1 

8x10 cm , It is in good agreement with the value of a re- 
ported by Ghoshtagore for oxide-silicon interface [20]. It is 
to be expected because the inert ambient as defined in Sec. 3. 2 
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contains a small amount of oxygen and therefore a thin PSG 
layer is expected to be present. The impurity therefore 
enters the silicon from PSG rather than from ambient. The 

value of ’a’ determined thus, should be expected to correspond 
to PSG-Si interface. 

The parameters A and D 2 were determined by using eqns. 
(2.16) and (2.21) for diffusion in inert and oxidising/ 
chlorooxidis ing ambients respectively. These equations have 
two unknowns, A and 02- The values of these parameters were 
determined by simultaneously solving the eqn. (2.i6)(or eqn. 
(2.21)) at two points on the experimentally determined impurity 
distribution (soe Appendix V). From the values of A and D 2 , 
thus determined, the theoretical impurity profile was generated 
using eqn. (2-16) (or eqn. (2.2l)). The theoretical profiles 
along with their experimental counterparts are plotted in Figs. 
3. 7-3. 24a. In most cases the two profiles are in good agreement. 
The details of tho experimental results are discussed in the 
following chapter. 
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CHAPTER IV 

EXPERIMENTAL RESULTS AND DISCUSSION 

4.1 PHOSPHORUS DIFFUSION IN PSG 

It was mentioned in Sec, 3.5 that in case^of strong 
oxidising ambient, the dopant concentration at silicon surface 
is time invariant and therefore the impurity distribution in 
silicon is given by, 


C«(z,t) == A.erfc ( . St . TXi ik ) 
2 2VD2t 

where, 

m Cq, 

A « 

1+ fit de* erf0(6) F(^^ ) 

2irD2 


(4.1) 


(4.1a) 



The values of were determined as described in 

Sec. 3.5. From eqn. (4.1a), it is seen that the value of 6 , 
and hence that of Dj^, can be calculated if and are 

known. was calculated with ni=10 [1], measured values of 

C ,B and estimated values of D, at each temperature. 

O' ^ 
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The ostimatod values of for OA samples (see Table 3.1a, 
Chapter III) are plotted against 1 oVt(°K) in Fig, 4.1. The 

least square fit to the exporimontal points is given by, 

_ 0-89 ^ 

^1 * 4.62x10 ( 4 . 2 ) 

rhi.“ value of activation energy of diffusion in tho above 
equation, l.c., 0.89 oV, is considerably lower than the reported 
values of 1.4~I,6 oV in literature [2-5], This is essentially 
duo to the difference in the nature of the experiments. In the 
oarliei reported studios, tho experiment consisted of diffusing 
phosphorus from n doped oxido source [3] or from a POCI 3 source 
[ 2 , 4 , 5 ] into thermally grown Si02. Tho diffusing phosphorus, 
presumably in form of [2-5], converted Si 02 into PSG and 

was estimated either from the movement of PSG-Si 02 interface 
[ 2 , 4,53 or from radiochemical profiling of diffusant in Si 02 
[ 3 ]. As pointed out by Eldridge and Balk [4], the activation 
energies reported by Soh et al [ 2 ] and Allen et al [5] (same 
applies to ref. [3]), actually refer to tho overall process of 
PSG film formation rather than to diffusion process alone. In 
contrast, in the present work, the PSG of uniform composition 
is grown during the diffusion process and the estimated 
essentially refers to the diffusion of dopant species in PSG. 

4.2 SURFACE ODNCENTRATION CONTROL BY AMBIENT 

Tho inpurity concentration at the silicon surface for the 
caso of diffusion in strongly oxidising ambient is given by, 



84 


~ A« srf c ( ) ^4 4 ^ 

2^02 

where A is given by eqn» (4.1a). A perusal of eqn. ( 4 . 1 a) 
points to an interesting and useful method of controlling the 
surface concentration i.e. by altering the diffusion ambient 
componitlon. According to eqn. (4. la), the value of A can be 
v«jried by orders of magnitude simply by varying B, i.e. the 
oxidation rate for the same value of Cq,D^ i.e. constant dopant 
concent ra tion in diffusion ambient. The oxidation rate can be 
controlled either by controlling the partial pressure of oxy- 
gen in the ambient or by changing the dif fusion/oxidation 
temperature. For low surface concentration, large values of 6 
(high oxidation rate) and D 2 (high diffusion temperature) are 
desirable, according to oqn. (4.1a). 

Figure 4.2 gives a plot of A vs 10‘^/T°K, as estimated in 
Sec. 3.5, for OA samples detailed in Table 3.1. The ambient 
composition during all these diffusions was kept the same. As 
would be expected, the composition of PSG grown in all these 
cases was also observed to be approximately the same 
(4 molo^ P 2 O 5 in Si 02 ) as already mentioned in Sec. 3.3. In 
Fig. 4*2, A varies almost by two orders of magnitude in 
the temperature range 950°C-1250®C. The lowest surface concen- 
tration is obtained at highest temperature, in agreement with 
the theoretical prediction. 




’ CwrV. ':n'd'“S^S^;KoA*-sSmp>., 
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In practical situations encountered in device fabrica- 
tion such as low concentration diffusion for base of p-n-p 
transistor, thyristor and hyperabrup varactor etc., this 
method can easily be adopted in preference to usual low 
temperature deposition/high temperature drive-in cycle. The 
growth of PSG during diffusion has advantageous side effects 
such as, (1) it getters the metallic impurities from bulk 
silicon, and (ii) the presence of PSG improves the uniformity 
of phosphorus diffusion over the surface. 

4.3 PHOSPHORUS PROFILES DIFFUSED IN OXIDISING AMBIENT 

The con^^arison of theoretical and experimental impurity 
profiles in Figs. 3.7-3.24, brings out an interesting fact 
that the profiles diffused in inert ambient (Figs. 3.11-3.14) 
are in excellent agreement with theory, while the same diffused 
in oxidising/chioro-oxidising ambient show a quantitatively 
minor but qualitatively consistent deviation from the calcula- 
ted profiles. In the latter cases, a pilo-up near the surface, 
a hump in the mid-region and a consequent kink in between the 
two regions is observed in almost all the cases to varying 
extent. A similar discrepancy between theory and experiment 
(minus the near surface pile-*up and hence, the kink) can also 
be seen in the experimental profiles reported by Masetti et al. 
[7,8], although they ignored it. 

The presence of the above mentioned anamaloue features in 
the profiles diffused in oxidising ambient alone, indicates that 
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the phenomenon is oxidation rel^toH tu • , 

uAxuouxun rej-ated. The pile-up essentially 

implies that the phosphorus is incorporated into silicon at 
PSG~3l interface, at a rate larger than that at which it can 
diffuse away into silicon. This is also corroborated by an 
infinitely large value of cx, hence that of reaction rate con- 
stant of phosphorus at the interface (see Sec. 3.5). The hump 
appears to bo an indication of the presence of non-Fickian 
movement, in addition to the normal Fickian diffusion, of 
impurity atoms, as postulated by Dobson [9,10] to explain the 
observed hump in the high concentration (> io^° cm~^) phos- 
phorus profiles diffused in partially oxidising ambient. 

Dobson [lO] assumed that oxidation of silicon induces a 
'convective m.ttorial flux', which includes both the self- 
intorstitials generated during oxidation as well as impurities. 
This flux was Tscumed to be proportional to oxidation rate 
and impurity concentration. The inclusion of this 'convective 
flux* in the model brings about an excellent agreement between 
theory and experiment as demonstrated by Dobson by deriving the 
equation for in^^urity distribution, assuming a linear PSG 
growth. For parabolic PSG growth also, as assumed in the model 
in Chapter II, the expression for impurity distribution in 
silicon can bo easily derived including the convective flux. 

The detailed analysis of the experimental results, on the 
basis of Dobson's model, however, has not been attempted 
because of the following reasons: (i) The physical basis of 
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the hypothesis is not clear. The mathematical representation 
of the 'convective flux' indicates that this flux is caused by 
a drivin., force [63]. This driving force is not clearly identi- 
fied. (il) The extent of deviation from theory, in this work, 
does not present a consistent picture i.e. the deviation does 
not show .1 systematic variation with temperature or crystal 
oricnUtlon. Therefore, the analysis of these results may 
not yield .any meaningful information. This inconsistency is 
boliijved to bo caused by two factors, (i) The smoothing of 
shoot rosistanco data tends to suppress the anamalous features. 
For dtfff’fcnt samples; the extent of data smoothing was 
diffutont, causing the suppression of anamalous features to 
varying extent, (ii) There is possibly some unidentified and 
non-reproducible surface effect contributing to this deviation. 
More oxperim».nts, preferably using direct methods such as SIMS, 
radiochemical analysis for impurity distribution evaluation, 
which may be free from data smoothing, are therefore required 
to be carried out to study those effects. 

It may be noted that in most cases, the quantitative 
agreement with theory is reasonably good. For the purpose of 
this work, the effective diffusion coefficient using the model 
of Chapter II has been calculated as described in Sec. 3.5, 

4.4 OED OF PHOSPHORUS IN SILICX)N-EXPERIMENTAL 

4.4.1 Diffusion in Inert Ambient 

The logarithm of the estimated diffusion coefficients are 
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ploucd against ioVT(°K) in Fig. 4.3. As can be seen, a 

straifjht lino can bo fitted to the experimental data. 

- eY 

= 1.74 e (4.5) 

Tho activation energy of diffusion, 3.56 eV, is in close 
agroomont with the other literature reported values 3.66 eV [ll] 
and 3.5 oV [7-8] respectively. The former was obtained by 
Fair [ll] by fitting a straight line to the (In D) vs. l/T 
oxperimontal data reported by a number of authors [12-15] while 
the later was reported by Masetti et al. [7-8], on the basis 
of tholr experiments on phosphorus redistribution during drive- 
in in Inert ambient. At 1250°C, tho diffusion coefficient is 
scon to be exactly same for both <100> and <111> crystal orien- 
tation in agroomont with tho observation of other workers, that 
the phosphorus diffusion is independent of crystal orientation, 
in inert ambient [ll]. 

4.4.2 Diffusion in Strongly Oxidising Ambient 

Tho In D vs IoVt^^CK) plots are given in Figs. 4.3-4. 5. 

The main features of phosphorus diffusion in oxidising ambient 
arc summarized below. 

(i) Tho phosphorus diffusion in oxidising ambient is enhanced 
in comparison to that in inert ambientj enhancement increasing 
with decreasing temperature. Tho diffusion enhancement ratio, 
increases from 1.3 at 1250^0 to about 4.7 at 950°C. 



Oif fusion Co-«ff. l€i»r/:i>cc- 



FlO.4.3 OKidflitton Enhanced Diffusion in <ioo>Silie 
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Fi6.4.4 Oriantatlon Dependence of Oxidation 
Enhanced Diffusion . 





Also plotted in Fig. 4.3, are the OED data of Antoniadis et al. 

[16] , obtained from redistribution of phosphorus during oxida- 
tion in dry O 2 * These are in good agreement with the results 
of this work. 

(ii) Tlic* phosphorus diffusion coefficient in oxidising 
ambient and therefore OED, is independent of crystal orienta- 
tion {see Fig# 4»4), This is in contrast to the behaviour of 
boron which exhibits an orientation dependent OED in silicon 

[17] . 

(iii) Tho OED is almost independent of oxidation/diffusion time 
( soc Fig* 4.b). 

In the following section the theoretical models proposed 
in litoraturo to explain tho OED as well as other oxidation 
related phert ^mcna, aro described and critically examined. 

4.5 OED OF PHOSPHORUS IN SILICON-THEORETICAL 

To understand and model tho phenomenon of OED, it is 
essential t-^ know tho identity and kinetics of the point 
defects Involved in OED as also the relationship between the 
point defect concentration and diffusion coefficient. 

4*5.1 Point Defects Involved in OED 

Since the diffusion in silicon proceeds via the motion 
of lattice point defects [18], the OED implies that oxidation 
introduces excess point defects in silicon* These point 
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defects, however, cannot be identified merely from the OED 
phenomenon, because the diffusion mechanism of phosphorus in 
silicon is still not known uniquely [18-22]. The possible 
diffusion mechanisms of phosphorus are considered to be 
interstltialcy [10,20], vacancy [ll] and a combination of 
both [ 18,20,21]. On the basis of OED alone, therefore, the 
point defects generated during oxidation could be either 
interstitials or vacancies. 

The oxidation is also found to cause the growth of 
stacking faults in silicon. The oxidation induced stacking 
faults (OSF* s) have boon determined, by TEM analysis, to be 
two dimensional extrinsic defects lying in <111> planes and 
bounded by Frank partial dislocations [23,24]. These extrin- 
sic OSF* s can grow either via absorption of self interstitials 
or emission of vacancies. The former implies the presence of 
supersaturation of solf-intorstitials and the latter that of an 
undersaturation of vacancies. Since both the OED and the OSF 
growth are observed to occur during oxidation under same condi- 
tions the point defects generated during oxidation are uniquely 
identified to bo self-interstitials [21,25]. 

4.5.2 Self Interstitials and Phosphorus Diffusion 

Since the oxidation causes injection of self-interstitials, 
the diffusion in oxidising ambient may be assumed to proceed 
either via interstltialcy or via a combination of both, the 
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interstitialcy and vacancy mechanisms* The diffusion coeffi- 
cient of phosphorus then may be expressed by a general eqn.[26], 


D 


ox 


D, 


21 


(-1 ) + Q 

ceq 

I 


2V 


(JL) 


(4*6) 


whore and are diffusion coefficients for interstitialcy 
and vacancy mechanisms, Cj, Cy the concentrations of self inter- 
stitials and vacancies during oxidation' and are the 

thermal equilibrium concentrations of self interstitials and 
vacancies respectively. 


As postulated by Hu [27] and experimentally verified by 
Antoniadis and Moskowitz [26], the recombination of vacancies 
and interstitials is a very slow process, limited by an energy 
barrier of about 1*4 oV. Under oxidising conditions, therefore, 
Cy sa Cy*^ and eqn. (4.6) can be written as, 
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°ox ° °2I °'- 

^ C, - C?'' 

* Dj + DriT 
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D, + D, 


(4*7) 


whero D* ■ Dot + is the diffusion coefficient in inert 
ambient and D^ denotes the oxidation enhanced diffusion. Since 
Cj is oxidation rat© and hence time dependent (this will be 
further discussed in following subsection), and are also 
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functions of time. It has been shown by Lin et al« [28], that 
when the diffusion coefficient is a function of time alone, 
the measured diffusion coefficient is simply the time-averaged 
value of actual time-dependent diffusion coefficient. In terms 
of time-averaged quantities, eqn. (4.7) can be written as 
[ 20 , 21 ] , 

“ Di ^ ^ j" (Ci-C®‘’)dt (4.8) 

“-I ° 

where measured value of diffusion coefficiont in 

the oxidising ambient. Eqn. (4,8) relates the diffusion co- 
efficient with the concentration of excess self-interstitials 
anti hwnce with the oxidation process. It may be noted that the 
depth dependence of Cj- has been ignored in deriving eqn. (4.8). 
This is a reasonable approximation to describe OED because jun- 
ction depths arc usually much smaller than the diffusion length 
of self-intors titiais in silicon. 

4.5.3 Kinetics of Solf-Interstitials Generated During 
Oxidation 

The formulation of the kinetics of the excess self- 
interstitials v/ouid require the knowledge of the exact relation- 
ship between oxidation process and the generation rate of self- 
interstitials and also the mechanisms of their motion and anni- 
hilation in silicon. Furthermore, the model should be able to 
account for tho other oxidation related ■phenomenon which is 
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obcorved under the same conditions i.e. OSF growth, both 
quantlt.itlvoly nr, well as qualitatively. Besides the experi- 
mental results detailed in Sec. 4.4, the other relevant experi- 
mental observations to be accounted for, by the model ,are the 

following . 

(i) The OSF’ s grow with time as (n = 0.77—0.84) while OED 

(of boron) varies as t ^ Since OSF growth rate ( d{lQgp)/dt) 

is proportional to excess self— interstitial concentration, both 

these observations are consistent with each other and point 
— 1/4 

towards t dependence of excess self— interstitial concentra- 
tion 

(li) Both the OSF growth rate and OED (of boron) are crystal 
orientation dependent and are larger on <100> than on <111> sur- 
faces [25]. 

(lii) Both OSF growth and OED increase with oxidation rate [25] 

(iv) The concentration of solf-interstitials at any point in 
silicon is diffusion limited. This is evident from the obser- 
vations that OED has boon reported to be a function of junction 
depth [29] and the delay in the onset of the growth of an OSF 
has been observed to be a function of the distance of the OSF 
from the oxidising surface [30]. 

Various models have been proposed to describe the kinetics 
of excess seif-interstitials [33-37]. The models proposed by 
Murarka [33] and Leroy [34] have already been critically 
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examined by Hu [25] and will not be considered here. The 
models proposed by Lin et al. [35], Tan and Gosele [36] and 
Hu [3/] are described below and examined in context with the 
various experimental observations. 

(a) The Model of Lin ct al. [35] 

An schematic representing the kinetics of excess self- 
interstitials in the Lin’s model [35] is shown in Fig. 4.6(a). 

It was pointed out, following Tiller [31-32], that the thermal 
oxidation reguires the creation of free volume at the inter- 
face at the rate of 1.25 , where is the interface velocity. 

This free volume may be supplied either by a flux of self- 
inti.rstitials leaving the interface or a viscous flow of the 
oxide. It was argued that if about 45j< Si at the interface is 
oxidised and the rest (55/.) is placed at interstitial sites in 
the oxide, the free volume will not be immediately needed. 
However, since oxidation is experimentally determined to be 
mure than 99^ complete, these excess interstitials placed in 
oxide will have to bo oxidised in other reactions. Based on 
the above arguments Lin ot al proposed a two step oxidation 
process* In the first step» the lattice silicon atoms at the 
surface undergo a reaction with oxygen to form Si 02 » 

Sij^ + ( 1-'Y)02 (l-Y)Si02 + YSij; (4.9) 

whore Si|^ and Si| represent lattice and interstitial silicon at 
interface respectively, Y is the fraction of excess silicon at 
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Interface and has a value between 0-0.55. A small part of these 
excess silicon atoms flow into silicon (j^ in Fig. 4.6(a)) and 
remaining excess silicon atoms form the flux into SiO^. The 
geni.'iation rate of excess self— interstials, as per eqn* (4.9) 
is related to JjL»T2 interface velocity V^, by, 

G = -gi = _ (^) = Jj_ + (4.10) 

whet'C r to SIO 2 to Si volume ratio: and o is volume of a 
silicon atom. 

The loss of solf-intorstitials via recombination with 
vacancies and jhsorption by stacking faults and other defect 
structures is approximated by an infinite sink of interstitials 
at a distance L- from interface. This distance is of the 
order of the r «ngo of OED (25pm as reported by Taneguchi [29]). 
Ignoring the convective flux duo to tho interface movement, the 
fxprussion for Is written as, 

h = 

s 

where Dj is diffusion coefficient of solf-intorstitials in 
bulk Si. 

In the second step, it is hypothesized that interstitial Si 
incorporoted into Si 02 in tho first step, are oxidised by the 
following reaction, 

Si + I O 2 i: SiO^ 


(4.12) 


100 


where x is an unknown parameter and possibly represents a 
mixture of several different SiO^^ type reactions. This SiO^ 
is assumed to be further oxidised to Si 02 later. 

The flux J 2 is determined by reaction (4.12) and is given 
by, 

•^2“ (4 - 4*) 

“ (ci - 4*) 

= 4 (4,13) 


Neglecting in ogn. (4.10) and combining with (4.13) 
leads to, 


(cj-q*) = 


r 

Q.Ki 


dXo 1 - f 


(4,14) 


It is further assumed that self “interstitial concentrations 
in SiOo and Si at the interface, are related with each other by 

(Cl 

an equilibrium stjgregation constant mj, 


,eq 




q-q^ 


(4.15) 


Using eqns. (4.14) and (4.15), the eqn. for excess self- 
interstitial in silicon, at the interface is given as, 
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On the basis of experimental results, Lin et al. proposed 

that the value of x is 1.2. For parabolic growth of oxide, 

X = 1.2 leads to the OSF growth as t*^*^ and OED variation with 
0 2 

time as t * , which compare well with their respective experi- 
mental values of (0.77-0.84) and 0.25 respectively. It is fur- 
ther speculated that ra^ in eqn. (4.15) is orientation dependent 
and ^ ''”'I<lll> * analogy with the experimentally 

observed orientation dependence of boron segregation in Si 02 ~Si 
system [l7,38]. The orientation dependence of OED and OSF 
growth rate was attributed to the orientation dependence of m^. 


Although this model explains some of the experimental re- 
sults mentioned earlier, it is unsatisfactory in some respects 
as described below. 


(i) An orientation dependent m^ would predict orientation 
dependent OED irrespective of the impurity type. While boron 
exhibits orientation dependent OED [17,39-41], the phosphorus 
OED in this work, has been observed to be orientation indepen- 
dent even under intrinsic diffusion conditions. 

(ii) It is implicitely assumed that the diffusion coefficient 
Dj. of self-interstitials is very large as compared to the 

range L of the self-interstitials, while correlating the 
concentration of self-interstitials at the interface (as given 
by eqn. (4.16)) with the time dependence of OSF growth rate and 
OED. The finite range of interstitials is attributed to the 
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loss mechanisms such as recombination with vacancies and 
absorption at defect structures and is essentially time 
invariant. The later experiments of Taniguchi [30] clearly 
demonstrated that the concentration of self-interstitials in 
bulk silicon is diffusion controlled. 

(iii) In the first step of the oxidation it was postulated 
that a fraction of interface silicon atoms become interstitial 
in Si02 via reaction given by eqn. (4.9), thus obviating the 
need of the free volume for the formation of Si02. However, 
when these are oxidised subsequently, the free volume is still 
required. It was pointed out by Tan and Gosele [36] that the 
model does not explain clearly, how this free volume require- 
ment is mot. 

(B) The Model of Tan and Gosele [36] 

Tan and Gosele [36] proposed a macroscopic rather than 
atomistic model for relationship between generation of excess 
solf-intorstitials during thermal oxidation. In this model 
the free volume requirement of the growing Si 02 is met by the 
viscoelastic flow of Si02 near the interface. The complete 
process is depicted in Fig. 4.6(b). It was argued that the 
growth of now Si 02 layer at interface amounts to insertion of 
oxygen atoms in bond centered position of Si atoms which would 
lead to a linear expansion of 30A in x,y,z directions. In z- 
diroction the expansion is accomplished easily as Si 02 surface 



104 


is free. In x~y plane, however, the Si02 layer is subjected 
to an elastic strain and associated compressive stress • t' . 
This stress causes the viscoelastic flow of excess material 
into Si02 (Fig* 4.6b(iii) and injection of self-interstitials 
into silicon simultaneously. The viscoelastic flow of Si02» 
is stopped as soon as the planar strain in x-y plane is fully 
released (Fig. 4,6b(iii)), The viscous flow rate is equal to 
the required volume expansion rate (and therefore oxidation 
rate) , The general relationship between viscous flow rate and 
stress T is given by, 

. dX 

Flow rate CO T CO (^) (4,17) 


where ' j’ is a function of viscosity of the newly grown Si 02 
near interface. Super saturation ratio of self-interstitials » ^rj 
is taken to be directly proportional to compressive stress ’r*, 
so that 
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r“ 




or 

(Cj. 




^^ 0 ) 1/3 
dt ^ 


(4.18) 


Tan and Gosole argued that ' j’ may take any value greater than 
or equal to unity. For well annealed amorphous Si02 J — 1. In 
contrast, the newly grown Si 02 layer has a lower viscosity so 
that j. > 1. It can be seen that for j = 2, (Cj-Cj^) » t ^ , 
which is the desired relationship to explain the time variation 
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of OoF growth and OED. It was argued that the scatter in the 
reported values of n (where OSF length oot*^) between 0.77-0.84 
is due to the variation in ’j' which is a strong function of 
impurity incorporated in the Si 02 near the interface. 

It is noted that the model ignores the diffusional effects 
oi self —inter stitials in bulk silicon and therefore cannot 
account for the observations of Taniguch [30] mentioned ear- 
lier. Furthermore, the question of orientation dependence of 
OSF growth rate and OED has not been addressed at all. 

The Model of Hu [37] 

This model is based on the following postulates, 

(i) The oxidation of silicon at Si-Si02 interface is usually 
incomplete, thus creating self-interstitials. The generation 
rate of solf-interstitials at the interface is directly pro- 
portional to the oxidation rate. 

(ii) At the silicon surface, there exists a certain density 
of active sites, possibly kinks of surface steps, which are 
capable of capturing the excess self-interstitials. The 
capture of self-intorstitials by these kink sites, thus 
causes surface rogrowth. On the basis of the calculations made 
by Burton et al. [42] which showed that the surface configura- 
tion entropy (and hence the density of the kinks) is a conti- 
nuous function of temperature, Hu [21] proposed that the ther- 
modynamic fluctuations can create the surface kink sites 
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at any temperature. He further explained that the density of 
kink sites is a function of crystal orientation and increases 
in the order <100>, <110> and <111>. Fig. 4.6(c) illustrates 
the process of capture of self-interstitials at kink sites. 

At site A the surface structure is restored by capture of two 
self-interstitials. 


(iii) The diffusivity of self-interstitials is finite and the 
loss of self-interstitials in bulk silicon due to recombination 
with vacancies and absorption by OSF’ s is very small. 


The distribution of self-interstitials in silicon, is given 
by the continuity equation. 


_ n 

FT “ 



(4.19) 


The boundary and initial conditions are. 


Cj(«.,t) == Cj(x,o) = 


(4.20) 


and 


[-D, 


dCi 




■cf)],, 


=0 


= Generation rate = G(t) 




(4.21) 


The solution to eqns. (4,19) - (4,21) is given by. 
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AC, 




TT 1 /O ^ +Vt 

^ ® e ^ erfc(x„+V-tJ 


n ' n' 
(4.22) 


where 


n 


n 


2VDj.t 
k? 


D, 


kj = Annihilation velocity (cm S“^) of self- 
interstitials at silicon surface 

Cj = Concentration of excess self-interstitials 

The growth rate of OSF was taken to be proportional to 
the < , the average excess interstitial concentration 

over the diffusion length, L^. The expression for < 
is derived from (4.22) and is given by. 


I o 


— — [e”^ erfc(l+ft„) 

X i-/-r rMfjL I* 


2ft 


n 


e ^ Grfc(ftn) + erfc(l)] 


(4.23) 


Hu demonstrated that the plots of ln( AC^.) vs* In(t^) 


and that of ln(<ACy>T ) vs. In(t^) have the slopes, of about 
(-0*25) and (-0.21) respectively in the range = 0. 1-1.0 
and of about (-0*4) and (-0,35) respectively in the range of 
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= 1 to 10. It follows that OED and OSF growth rate vary as 
t and t respectively (in the range t^ = 0.1 to 1), 

in agreement with the experimental results. Hu argued that 
most of the experimental results reported cover only about 
a decade of time span and are likely to correspond to range of 
= 0. 1-1.0. The small variation (n = 0.77-0.84) observed in 
the time dependence of OSF growth rate is merely due to the 
different value of t^^ in the various experiments. The only 
variable in eqn. (4.22) that depends upon crystal orientation 
is kj. and is responsible for orientation dependence of OED and 
OSF growth rate. Since the <111> surface has a larger density 
of kink sites than <100> surface, k-j- would be larger for <111> 
surfaces, resulting in smaller OSF growth rate and OED for 
<111> surfaces, as experimentally observed. The assumption of 
the diffusion controlled motion of self-interstitials is also 
consistent with the observations of Taniguchi [29,30] mentioned 
earlier. This model, thus, is qualitatively and to some extent 
quantitatively consistent with, the experimental observations 
reported earlier. It will be shown below that this model can 
also explain the orientation and time invariance of OED observed 
in this work. 

In the model described above, k^ is considered to be 
dependent on the density of kink sites. It is proposed that 
kj would also depend upon the impurity present at interface 
through its segregation coefficient in Si-Si02 system. Since 
the segregation coefficient of phosphorus in Si-Si 02 system is 
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large (m = 10), it may be assumed that phosphorus will tend to 
pile~up on silicon side of interface, and effectively compete 
with s elf "in ters ti ti als in surface regrowth process. It 
is proposed that phosphorus is captured at kink sites in 
preference to self-interstitials. Since density of kink sites 
is estimated to be k lO"-^ cm”^ [21], a phosphorus concentration 
of about 10 cm at the interface would be sufficient to eli- 
minate the capture of self interstitials at the interface. This 
moans that k^ ts o in eqn. (4.22). This situation will load to 
the following effects. 

(i) Cj.(x,t) and hence OED, would be orientation as well as 
diffusion time independent, for case of small junction depths 
such as Xj << 2VDjt, 

(ii) Larger OED than expected for the case of kj > o would be 

observed since CjCx,!) in eqn. (4.22) reaches its maximum 
value as ** o* the present work, since the surface 

j[0 —3 

concentration in almost all cases is of the order of 10 cm , 
all those effects may bo expected to be present. 

It follows from the above argument that for group V impuri- 
ties (P,As,Sb), which have m >> 1, the concentrations of the 

1 O 

order of 10 cm and above will lead to orientation and time 

16 *3 

invariant OED. For lower concentrations, say <10 cm , the 
OED (as well as OSF growth rate) should be both orientation 
and time dependent. Unfortunately no systematic steady has been 
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reported on orientation dependence of OED of P,As and Sb. How- 
ever* it has been reported by Antoniadis and Moskowitz [26] 
that for concentrations of the order of 10^^ cm"^, P and As 
indeed exhibit a time dependent OED. Furthermore, for boron, it 
would be expected that OED is orientation dependent even for 
comparatively higher boron concentrations at the interface, be- 
cause boron tends to segregate into Si 02 , thus causing the 
surface regrowth to occur dominantly by capture of self-inter- 
stitials* In literature, the orientation dependence of boron 
has been reported for impurity concentrations ranging from 
10^”^-10^^ cm'“^ [17,39-41]. 

Conclus ion 

The review of the models, proposed to describe the kinetics 
of excess self-interstitials leads to the following conclusions* 
(i) Only the model proposed by Hu [37] is consistent with all 
the experimental observations on 03F growth and impurity OED in 
silicon (ii) The surface regrowth process in Hu’s model is a 
function of the segregation coefficient and concentration of the 
impurity present at the Si02'-Si interface. It may, however, be 
pointed out that in all these models the Si 02 --Si interface has 
been assumed to be a source of self-interstitials which is con- 
sistent with the phenomena of OSF growth and OED. At high tem- 
peratures (>1150*^0) and long oxidation times, however, OSF shrin 
kage and ORD of phosphorus are observed [43,44] which point to 
the fact that under these conditions the interface acts as a 
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sink for self-interstitials [43] . This situation has not 
been considered in those models. The reason probably is that 
the precise mechanism of the oxidation process and generation 
of excess self-interstitials is still not well understood. A 
comprehensive model to describe the generation kinetics of 
excess self— interstitials at the interface is, therefore, yet 
to be developed. While the two basic features of Hu’s model 
i.o. surface regrowth process and finite diffusivity of self- 
interstitials in bulk silicon have to be retained in any model, 
the generation process of self-interstitials at interface needs 
to bo reformulated. However, as shown in the foregoing discu- 
ssion, Hu’s model in its present form satisfactorily accounts 
for all experimental observations on OSF growth and OED and 
hence can be considered to be valid under corresponding oxidi- 
sing conditions. 

4.5.4 Modelling of OED 


Following Hu’s model for kinetics of excess self-intersti- 
tials, generated during oxidation, the equation for oxidation 
enhanced diffusivity factor ’G’, for the present work (i.e. 


kj. = o) can be written as (using eqns. (4.8) and (4.22)), 


<D > 

r - ^ 1 

Q - D] 

i 


Dj^ — A-f> 

^21 I 




(4.24) 


It is evident that a plot of In G vs 10 VT(°K) would be 
straight line. In the plot of Fig, 4,7, the experimental points 




Fig. 47 Oxldtttfon Enhancad Olffusivily Fader 
<100> Silicon. 
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are best described by a straight line, given by the equation, 


G = 1.246x10 


-7 


1.84 
J kT 


eV 


(4.25) 


Equation (4.25) relates the oxidation enhanced diffusi- 
vity factor with temperature. This equation is valid only 
for diffusion in 90/. 0^ ambient. 

In eqn. (4.24), G is a function of oxidation conditions 
only through Aj. By evaluating G for different oxidation 
rates, it is possible to check whether or not the interstitial 
generation at interface is proportional to oxidation rate as 
assumed in the model of Hu [37]. 

The value of parameter at low phosphorus concentra- 

tions (^ 10^^ cm ^) can be calculated, if the oxidation/dif fu- 
sion time dependence of G is known. The only data available 
on oxidation time dependence of phosphorus diffusion coeffi- 
cient are those reported by Antoniadis and Markowitz[26]. 

Their data, however, are reported only at one temperature i.e. 
1000°C. More experimental data are required before a detailed 
evaluation of this parameter is attempted. 

4.6 THE EFFECT OF CHLORINE ON PHOSPHORUS DIFFUSION 
4.6.1 Experimental Results 

The phosphorus diffusion coefficient in <111> and <100> 
silicon is plotted against lotT(°K) in Fig. 4.8 and Fig. 4.9 
respectively. It is seen that 
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(i) The presence of chlorine, (produced by decomposition of 
TCE [ 45 ]) in the ambient causes a reduction in the diffusion 
coefficient at all the temperatures. 

(ii) In <100> silicon the reduction is less than the same in 
<111> silicon. 

(iii) The effect of chlorine on the diffusion coefficient 
decreases with decreasing temperature. 

(iv) At high temperatures (> 1150°C) the diffusion coefficient 
in chlorooxidising ambient is lower than its value in inert 
ambient. In <100> silicon, the diffusion coefficient at 1250^0 
is reduced by lO/. and in <111> silicon by 35?«i, with reference 
to its corresponding value in inert ambient. 

In order to explain the above results it is essential to 
know the distribution of chlorine in Si 02 --Si system during 
chlorine oxidation and also the nature of interaction of 
chlorine atoms with the point defects involved in impurity 
diffusion. 

4.6.2 Chlorine Distribution in the SiO^-Si System 

Recently, there has been considerable effort invested in 
understanding the effect of chlorine on thermal oxidation 
because chlorine in thermally grown oxides (i) neutralizes the 
positive alkali ions in the oxide [46], (ii) provides better 
reproducibility in breakdown voltages [47], (iii) and reduces 
stacking faults [48]. A brief review of the status of chlorine 
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FIG.4.8 The Effect of Chlorine on Phosphorus 
Diffusion In <ni> Silicon. 
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FI6.4.9 The Effect of Chlorine on Phosphorus 
Diffusion in <ioo> Silicon. 
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oxidation in the fabrication process of LSI's is given by 
Hattori [49]. 

Rather a few investigations have been reported on the 
kinetics of chlorine during the thermal oxidation of silicon 
in chloro-oxidising ambient [50-53]. The chlorine distribu- 
tion in Si 02 in these studies was experimentally determined 
using nuclear backscattering [50], SIMS [51] and Auger sputter 
prof iling/SIMS [52-53]. The main observations are as follows : 

(i) The chlorine atoms segregate near the interface in Si02» 

where these are probably incorporated in various Si-Cl-O-H 

complexes. The peak in the chlorine concentration profile 

occurs very near the interface on Si02 side. There is a 

region of 100-200°A thickness near the interface where chlorine 
is accumulated in high concentrations [50-53]. The concentratii 
of chlorine in the bulk Si 02 is about two orders of magnitude 
lower than the peak concentration near the interface. 

(ii) The chlorine concentration at interface is a strong fun- 

13 

ction of tenporature. The chlorine level varies from 6x10 cr 
for 900^0 oxides to about 1.7x10^^ cm~^ for 1100°C oxides ; 
[50-53]. The clorine concentration also is a weak function of 
oxidation time [51,52] and increases linearly with chlorine ; 
partial pressure in ambient [5l]. ! 

(iv) For identical oxidation conditions, the level of chlorim 
incorporated in Si 02 on <111> silicon is larger than in Si 02 o| 
<100> silicon. This experiment, however, was carried out only | 
for oxides grown at llOO^C [51] and 1000 C and 1100 C [52] . | 
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( v) The linear rate constant of oxidation is increased in 
presence of chlorine [54-56], 

4,6,3 Interaction of Chlorine with Self-Interstitials at the 
Interface 

It is now well established that oxidation in chloro- 
oxidising ambient causes the retarded growth or elimination 
of OSF [57-59], The reduction in OED of phosphorus and boron, 
in chloro-oxidising ambient, has also been reported by Nabeta 
et al, [60], As pointed out in previous section, both these 
phenomena i,e. OSF growth and OED are caused by the excess- 
self-interstitials generated at the interface during thermal 
oxidation process. It thus follows that the chlorine 
incorporated at the interface during thermal oxidation causes 
the interface to act as a sink for self-interstitials, 

A two step mechanism has been proposed by Murarka [61-62] 
that possibly reduces the concentration of self-interstitials 
near the interface. In the first step, it is postulated that 
the interface chlorine atom reacts with the lattice silicon 
atom to form SiClu, by the reaction, 

SIl + 5 CI2 Z SiClu + Siy (4.26) 

The silicon vacancy generated in this reaction recombine 
with the self-interstitial present at the interface, 

Sij + Siy it Si^ (4,27) 
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In eqn. (4.26), the interface chlorine atom has been 
assumed to react with lattice silicon atom alone. While it 
is consistent with the observed increase in linear rate con- 
stant during chlorine oxidation [53-55,62], it is unlikely 
that the free self-interstitials generated at the interface 
would not directly react with chlorine atoms. The exact 
chemistry of the interaction of chlorine atoms with various 
species such as self-interstitials, 0 and H-atoms present 
at the interface is not yet known. It appears, however, rea- 
sonable to assume that in general the interface chlorine 
atoms react with self-interstitials forming various immobile 
SiCl^ or Si-Cl~0-H complexes [52]. The Si02“Si interface 
containing chlorine, thus, acts as a sink for self-intersti- 
tials generated at the interface. 

4.6.4 Discussion 

In general, the loss rate, R, of self-interstitials at 
the interface duo to reaction with chlorine atoms may be 
expressed as, 

(4.28) 

where k^^^ is a constant, and are the concentrations of 
chlorine and excess self— interstitials at the interface res- 
pectively. Considering that Cqj^ is a weak function of time 
and can be treated essentially as a constant [52], it is easy 
to formulate the effect of chlorine on Cj-, in Hu‘s model 
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described in previous section. The effect of the presence of 
chlorine, as expressed by eqn. (4,28) can be included in the 
boundary condition given by eqn. (4,21), The modified eqn. 
(4.21) is, 

dC-r . 

(4.21a) 

The expression for excess self-interstitial distribution 
at silicon surface is determined to be (see Appendix VI), 

X/2. t* 

(Cj—Cj*^) I ~ ^Cj(o,t) ~ A j ( q“) e erfc(Y^t^) — 

t' 

[l-e ^ (4.29) 
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The experimental results are now easily explained on the 
basis of the qualitative discussion presented in subsection 
4. 6. 2-3 and its quantitative representation in eqn. (4.30). 

It is seen that the second term on RHS of eqn. (4.30) de- 
creases as Cq 2 increases and Aj decreases. The orientation 
dependence of the effect of chlorine on diffusion coefficient 
is essentially due to orientation dependence of Cq^. As 
mentioned in subsection 4.6.2, the is larger for <111> 
surface* This is consistent with larger reduction in diffu- 
sion coefficient for the case of <111> silicon, than for <100> 
silicon, in identical ambients. The decrease in the effect of 
chlorine on with decreasing temperature is due to the 

fact that % sharply decreases with decreasing temperature, 
as mentioned in Sec. 4.6.2. 

Lastly we come to the observation that at high temperaturei 
(> 1150°C) , the diffusion coefficient in chloro-oxidising, ambi- 
ent is lower than its value in inert ambient. This essentially 
implies that the second term in eqn. (4.30) is negative, signi-. 
fying a situation when at the interface. It is point- 

ed out that if is entirely due to the vacancy mechanism as j 
believed by some investigators [ll], the situation at | 

the interface (hence in the impurity diffusion region) will 
either leave Dj_ unchanged or will cause an enhancement in D^. ! 

The former situation would arise, if the process of recombina- : 
tion of self-interstitials and vacancy is slow [26], so that 
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vacancy concentration in the diffusion region, (hence re- 
mains nearly equal to its thermal equilibrium value. On the 
other hand if the recombination process is fast, the concentra- 
tion of vacancies in the diffusion zone will rise, causing an 
enhancement in D^. The observed reduction in diffusion coeffi- 
cient below D^, in chloro-oxidising ambient, therefore, uniquely 
proves that in inert ambient phosphorus diffuses via a pure/ 
partial 'interstitialcy mechanism. 

The fractional contribution of the interstitialcy mechanisn: 
(if a partial interstitialcy mechanism of phosphorus diffusion 
is assumed) to the phosphorus diffusion in inert ambient may 
also probably be experimentally determined by evaluating the 
effect of chlorine at higher chlorine concentrations. 


It can bo seen from eqn. (4.29) that for large such 
that Cci) » 1 , AC 3 . - which signifies that the 

g- 

interface is a perfect sink for self-interstitials. For such a 
situation the expression for Dox,cl reduces to. 


^ ^ox,cl^ 


D. 


D, 


21 


(4.31) 


If diffusion mechanism is partially interstitialcy, may 
be written as, ^ sufficiently high 

<D ,>wtll saturate to a minimum value, equal to The 

fractional contribution of interstitialcy mechanism, then can be 
calculated by the equation, 
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f = [1 _ (4.32) 

i 

It may be pointed out that in the present work, the volume 
percent of N2» bubbled through ICE (the source of chlorine in 
the present experiments) was merely by, of the total ambient 
gases. This can probably be increased by a factor of two 
without causing any damage to the surface [52]. This method 
may allow determination of f^, if f^ is less than, say 0.6. 



CHAPTER V 


surw^RY 


It is well established that the phosphorus diffusion in 
silicon is influenced by the diffusion ambient i.e. the 
presence of 0^ and Cl^ in the ambient. The objective of the 
present work is to understand and quantitatively model the 
ambient effects on the phosphorus diffusion. 

Several workers have systematically investigated the 
effect of oxygen on phosphorus diffusion. Their experiments 
consisted of evaluating the redistribution of deposited phos- 
phorus (either by CVD or ion implantation) in oxidising ambient. 
The analysis of these profiles for extraction of diffusion para- 
meters is rather conplicated, approximations are to be made for 
tho initial deposited profile and the finer details of the 
nature of experimental profiles are likely to be overlooked. In 
the present work, tho constant source diffusion using PH 3 , has 
been utilized to investigate the effect of 0 ^ and Cl^ on phos- 
phorus diffusion in silicon. 

A model for constant source diffusion in oxidising ambient 
is presented. Since in general, the surface concentration in 
silicon is time-variant, a reaction rate limited flux boundary 
condition has been assumed at the interface. It is shown that 
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in extreme cases such as negligible PSG growth and infinitely 
large reaction rate constant, the model yields Smit’s equation 
and Grove's equation respectively, for impurity distribution 
in silicon. For the case of small diffusion lengths, the model 
leads to a series solution. 

To investigate the ambient effects, the phosphorus diffu- 
sion was carried out in inort ambient, 90 ^2 

4.5 (N 2 +TCE) respectively, at various temperatures in the 
range 950°C'-1250°C. The sheet resistance profiles were mea- 
sured by Vander Pp.uw/ four— point probe and anodic stripping 
of silicon layers. The impurity distribution in silicon was 
determined by incremental sheet resistivity method. The 
diffusion coefficient of phosphorus was calculated by solving 
the appropriate theoretical equation for impurity distribution 
at suitable points on experimentally determined profile. The 
main results are the following. 

(i) For the case of strong oxidising ambient, the model pre- 
dicts a strong dependence of surface impurity concentration in 
silicon, on the ratios of PSG thickness (X^) to the diffusion 
lengths in PSG = 2 fD 3 ^t) and silicon (L 2 = res- 

pectively. The low surface concentration requirements are a 
large value of X^/L^ (i.e., strong oxidising ambient with low 
dopant partial pressure) and small value of (high tem- 

perature). It is verified experimentally also.' It was obser- 
ved that for the identical oxidising diffusion ambients, the 
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surface concentration varies by two orders of magnitude from 
its lowest value at 1250°C to its highest one at 950°C. 

(ii) It has boon observed that the detailed nature of experi- 
mental profiles, diffused in strong oxidising ambient deviates 
from theory, qualitatively. The experimental profiles show a 
pilo-up near the surface and a hump in the bulk. The pile-up 
is attributed to the infinitely large value of reaction rate 
constant of phosphorus at silicon surface, which causes the 
incorporation of phosphorus into silicon (at the interface) at 
n larger rate than that at which the phosphorus can diffuse 
into bulk silicon. The hump essentially indicates the presence 
of non-Fickian impurity motion, the origin of which is not 
known nt present. The deviation in most cases, however, is 
very small and has been ignored in calculating the diffusion 
couff icionts- 

(ill) The amount of OED observed in 90/, O 2 diffusions is 
approximately the same as reported by other investigators in 
dry O 2 ambient. The OED has been found to be independent of 
the crystal orientation and oxidation/dif fusion time. The 
models reported in literature, to explain the oxidation induced 
phenomena, have boon discussed in context with the experimental 
observations of this work coupled with the other related obser- 
vations reported in literature such as t ^ dependence of OSF 
growth rate and OED, their orientation dependence etc. It is 
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shown that only the model proposed by Hu can account for all 
the experimental observations reported so far. It is suggested 
on the basis of the observations of this work that the surface 
rogrowth process in Hu* s model depends upon the type of impurity 
as well as on its concentration at the interface. Finally, it 
is observed that more experimental data are reguired for quanti- 
tative evaluation of parameters in Hu*s model. 

(iv) The presence of chlorine in the diffusion ambient has been 
observed to cause a reduction in OED. The effect of chlorine is 
found to bo dependent on crystal orientation and temperature. 

The reduction in OED is consistent with the reported observation 
of retarded growth/elimination of OSF in chloro-oxidising ambient 
and is attributed to the fact that the chlorine segregates in 
Si02 at the interface and the interface containing chlorine acts 
as a sink for self-interstitials generated during oxidation. 

Hu’s model is extended to take into account the effect of chlo- 
rine on OED. It is shown that the effect of chlorine on OED 
varies inversely as chlorine concentration at interface, Cq^^. 

The observed orientation dependence of the effect of chlorine 
is explained to bo due to the reported orientation dependence of 
. Comparatively small effect of chlorine on OED at lower 
temperature (< lOSO'^C) is attributed to the reported sharp 
decrease in with decreasing temperature. The larger genera- 
tion of self-interstitials at these temperatures further obscu- 
res the effect of chlorine on OED. 
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The most important observation is the reduction of the 
phosphorus diffusion coefficient in chloro-oxidising ambient 
below its value in inert ambient. It leads to the conclusion 
that in inert ambient, the phosphorus diffusion in silicon 
proceeds via either partial or pure interstitialcy mechanism. 

Lastly it is suggested that these experiments may be 
repeated for higher chlorine concentration in the ambient to 
experimentally determine the fractional contribution of inter- 
stitialcy mechanism to phosphorus diffusion in silicon. 
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APPENDIX I 


Applying boundary condition given by eqn. (2.3), to the 
expression for Cj_(x,t), gives 

C^io.t) = C = Z |- n 

n=o 

The coefficients A^^ may be eliminated by comparing the 
coefficients of like powers of on both the sides 


Cq = erfc(o) 

— a + A_ 

0 0 

so that Aq =: (Co-a^) (2) 

For odd terms of the series of eqn. (1), 
o = a^.h^(o) 4- Aj^«i^ erfc(o) 

= A^.i"^ erfc(o) 

so that, ^ ® ^ ^ ~ **• 

And for even terms of the series, 

0 = ay^.hj^(o) + erfc(o) 

-105)1 • ®n ^ 2M^(| i"I) 

so that, * -2'^.nl.a^ ! n = 2,4, ... 
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Now since 


dx^(t) 

~dt ' 


■\rR T Di 

*CjL(x,t) = ~ 2ft* ^ •Cj_(x,t) 


the expression for the L.H.S. of impurity conservation condition 
(2.5), can be expressed as, 


*-°l ST" ^ dt '^lJx=X (t) = 27 t ao'^^V^o'- ‘2 6.erfc(6) 


TDi 


2e' 


+ CDit)"/2.3^ {2n.h^_^(6)+2 6 .h„(6) 


fit 


2 «» 

+ erfc(6) 


n=2 


26 .i*^ erf c( 6) }] 


3 i 0 e 9 


and 


^n+1^^^ 2x.h^(x) + 2n.h^_j_(x) 


2(n+l) erfc(x) = i*^”^ erfc(x) - 2xi’^ erfc(x) 


We hove, 




t°l TT FT* <^l^x=X„(t) = 5ft t26 .a^+CC^-a^). (2fi.erfc(6) 




2 x ^ 


erfcCa) ] 


(7) 
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APPENDIX II 


Transforming the dependent variable 02 ( 2 , t) in eqn. 
(2.2a), to C(z,t) = 02 ( 2 , t) + Cg, the resulting diffusion 
equation and boundary/initial conditions are, 


D 


d^c 

2^ 


4 * 


rfB dc dc 
2tt = -St 


( 1 ) 


C( 2 ,o) = C(oo,t) = 0 


( 2 ) 


C|f - = -m.a.CjL(X^(t),t) - aCg (3) 

The solution to eqns. (1) - (3) can be obtained by using 
Duhamel’ s theorem and is given as, 

i t 

C(z,t) ~ f {-aCg-m.a.C^CX^CX) , X) } .N(z, t- X)dX (4) 

whore N(z,t) is the solution to eqn. (1) with C(z,t) replaced 
by N(z,t) for the following boundary/initial conditions. 


N(z,o) == N(<»,t) == 0 

(2a) 


(3a) 


Now if a transformation of 1 = -^ -aN is carried out in 
eqn, (1), with C( 2 ,t) replaced by N(z,t), the resulting diffu- 
sion eqn* in ^ and boundary/initial conditions are, 
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D ® 4 . ifB d'i 

°2 ^ + 5tt H 

(5) 

f ( Z,o) = ^(oo,t) = 0 

(6) 

W(o,t) = 1 

(7) 

Solution to eqns. (5) - (?) is, 


S(2,t) LS = p.erfc(£talM) 

(8) 


erfc(£ g _ t _.,..) 


aiTDj 

where 

" = erfc(^l^) 

2fD2t 

Thus, 

|li ctN = ^ = p.erfc( — (9) 

2fD2t 

General solution to eqn. (9) is, 


N( z, t) 



az , az 
e + e 


} p erfc(i^±i3tai-).e-“^'.dr' 
0 2f D^t 


where Kji is an arbitrary constant. Since N -* o as z -*• » , it 
is evident that Kj^ - 0. Thus, 


N(r,t) = p.e“^ erfcC^^^-SCSi) .dr' 

0 ZfD^t 

= - S. . [erfc(^±E>CB±) - sd(2+i'V'Bt)+a2D2t_g^^^(Z«!£Bt 

“ 2TD,t ZVDgt 

^ + aTD,t)] 

g(z*t) 




oc 


( 10 ) 
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where , 


g(z,t) = + afD^t)] 


2VDot 


2fD2t 


( 11 ) 


The expression for C2(z,t) now can be written as, 


C2(z,t) = -Cg + / {Cg+m.CjL(XQ(X),X) }.g(z,t-X)dX 

Substituting for Cj_(X^, X ) , we have, 

*fc 

C2(z.t) = -Cg + ^ } [Cg+m tajj+(Cjj-aj,)erfcC 6) ) 

+ Z m d"/ 2 h (a) x"/2 

1 X r* n 

n=i 

« m 2’^.nl .a^.i”^ erfc(6) x'^/^] ,g( z, t-X)dX 

n=2 ^ 


“Cg+[CB+m. a^+CCQ-aQ) .erfc(6) ].g(z,t) 
,n/2 


+ E m.D‘j;'^.a^.h^(6). ^ X^^^.gC z, t-X)dX 
n*l 


E m.Dj^^.2’^.nl.a^.t’^ erfc(6). g( z, t-X)dX 

n«l o , , 

(12) 


Now since! 


^(z»l)l = S(t) = [-o£.8.e“^^'^^^®'*^^‘*'“^°2‘^.erfc(^i^^^ + 

¥z *z~o ^ ^ *■ 2fD2t 

aVDot)] 


2"'-'z=o 


-a 


2Y‘D2 



and 


A“7 


g(o,t) = p.erfc(£!^ + ccfD^t) 

^^2 a^D, 2 


1 + 


slii 


a 


the R.H.S. of impurity conservation condition (2.5) can be 
written as, 

^^2 d'z""" “*■ ift ^2 ^z=o ~ ■■ erfc(6)}] , 

[D 2 ^(t) + ^ (1 + ®-^)] 

OO y *t / 

+ S n, Df2 h„(s)[D2.|^ I |f U^o-dX 


n=l 


+ f X )dX ] 


- m 0^/2 2".nl.a„ l" erfc( ) lflz=o-^ 

4" S ^^^^•9(o,t~ X ) dX 

2ft ^ 0 

^ — [-rfB.Cg+[CQ+m CaQ+CCQ-a^) erfc(6)} ].[( 2 D 2 ft + 

2 "^ *t# 

El&).S(t) + r^B] 

t 


n=l 


I m a^ h^(6) .[(2D2ft+ ^ t-X)dX 

+ rfB t’^/^3 

m Dy^.2'^.nt.a^.i'^ erfc(6) [(2D2ft+ ^ S 


n~2 


S(t-X)dX + rfB.W^]] 


(13) 
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appendix III 


Using eqns. (2.12) and (2.14), the impurity conservation 
condition of eqn. (2.5) can be written as, 


26 126 .erfc(6)- 




2 ,oo 

- E 
n=2 




,n/2 


•2r.6 .Cg+[CQ+m {aQ+(C^-aQ)erfc( 6) } ] . [S( t) .Q( t)+2r. 6t] 

a } ^"/2.S(t-X)dX +2r.6.W2 


+ E (m.D" 
n=l 


— E ( .2 .nl . i erf c(6 ).a^)* Q(t). J . S( t— \ ) .dX 

n ^2 o 

+ 2r. 6 (1) 

where Q(t) - — (2D«V't + ^)= 2M.'V'Dot + 

VD^ / a 2 a 

pi = YD^/fDi 


and S(t) is as defined in Appendix II. 

The R.H.S. of eqn. (1) is rather complicated and it is 
difficult to determine the value of a^^' s for a general case. 
Some special cases, however, are considered below. 

Case I : 5 fS o 

For this case, the odd terms of the series on both the 
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sides become zero and even terms cancel each other on their 
respective sides. The L.H.S. then reduces to a constant 

Comparing the constants on both tho sides in eqn.(l), 
then leads to 


^o"”^o ^ 


(2) 


or 


^0 ~ *^o 


p(t) = mCj^CX^Ct) ,t) = mC, 


(2a) 


Case II : afD^t >> 1 


For afD^t >> 1, S(t) 


0 •§. 


•r^B/4D2 


and 


|t / x‘^/^.S(t-X)*dX = - - 

0 nfDr 


Furthermore, approximating 

af D^t 


^b|4D2 


r (2^) 

’ ‘rc^) ■ ^ 

0 in eqn. (l) and 


comparing the coefficients of like powers of *t’ on both sides 
of eqn. ( 1) , one gets, 

- 

2.6* a^+C CQ^a^) * { 26 • erf c( 6 ) “ ^ = 2r • 6 *^3 ^ Cg"^niL a^ 


+ (C^-a^).Grfc(6)]} . I2r6 - ^ B/^^2} 


W 


(3) 


= a„[m.D;/2.h^(5) { rfB - yf irDg.p.e'’^ 


r 

r(^) 


ill ^ ^ 1^3>&> • • • 


(4) 
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and, 

- 2'^.nl,i''^ erfc(6)}. i rfB -~ 

fn 


erfc(6)] 

TD^.p.e^ 


= Ch (6) 

1 n n' ' 

B/4D2 LlIII^ 


n =2,4,6, • . • 


(5) 


From eqn. (4) and (5) it is clear that a = 0 for all 

n 

2. From eqn. (3), the value of and hence the expression 
for aQ+CC^-a^) erfc(6) is derived to be, 


A = m[a^+(E^-a^) .erfc(6)] = 

mEc^Ce"’ * /Y'TC)-(ii.p/Y’ii).e”^ bAd^ (6).Cg] 

— ^ — ^ ( 6 ) 

( o”* ^ /^%)+ 6. erf (6)+m.r. 6. erf (6 ) .F(r,B,D 2 ) 


where , 

F(r,B,D2) = -1 + 1 

ftt(rVB/2VD2) .e^ ®/'^°2.erfc(rfB/2V'D2) 

C6a) 


Also, for “fDgt >> 1, 

g(z,t) = p,[erfc(^~^) 

2fD2t 


Y'n;,aY'D2t 


-] 


(7) 


p(t) =* mCj^CX^C t) ,t) = A 


(7a) 
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Case III : afV 

With this approximation, S(t), Q(t) are written as, 

S(t) B -a.p.e“^V.erfc(alfD2t); Q(t) » 2 urD 2 t and 

[s(t) Q(t)+2r6 ] fSS(t).Q(t) 


Ij- / K"^^.SCt-X)dX 




where 


A 9 r 1 S( t) 1 

Si(t) = -St S(t) = o. f 

Substituting eqn. (9) in eqn. (8), 

; / x"''2.S(t-X).dX + o;2^YD2. 


rc^) ^ 
r(^) 


CIO) 


X c/'+') nft^ and the definite integral of 
Substituting for S(t), QKV 

vie,-. VI V.> 

n/2 ni 

i" erfc(6).a„.lCt) (U) 
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where, 

I(t) = Q(t). |r/ ^ 

o 

- 9n„2 n3/2 ir+ An/a er. 2 ^ 

- 2p,a D2 .yt S X ' .S(t-x)dx+ 2 n,a^pD^. -5- .t ^ 

0 2 p(n^) 

- 2ti(xp y'D2.t ^ + 2r6 (12) 


The L.H.S. of eqn, (1) remains unchanged in the situation 
considered# Letting t -► o"^, on both sides of eqn. (1), with 
its RHS replaced by eqn. ( 11 ), one gets, 


= -•2r. 6 .Cr 


(13) 


26 aQ+(CQ-aQ).[26 .erfcCa) - 

From eqn. (13), first the expression for a^, thence the 


same for A as defined in eqn. (6), is derived as, 

- x2 

C (e /fit - r. 6 .erfc(6) C, 


m 


'B 


(14) 


(e" ^ /Vic) + 6.erf(6) 


Substituting eqn. (13), (14) in eqn. (1) with RHS replaced 
by eqn. (11), the eqn. (1) can be written as. 


00 


n 


/2 


. .t’^/2 « ^2 [D5|/^.2’^‘^^.(n+l)l. 

“ r,=0 


.n+1 

1 


n=2 

erfc(6) ]aj^t 
,n/2 


n/2 


2A.lifD«t* S(t) + 2 '^l ^ hj^(6) aj^.I(t) 

z n=l 

- m erfc(6 ) .a^.l( t) 

n=»2 


( 15 ) 
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The value of a^' s may now be deferrained by successively 
differentiating both the sides w.r.t. ft and letting t - o. 

It can be seen that the series terms on RHS of eqn. (15) will 
always tend to zero as t -* o after each differentiation. The 
coefficients s are given by, 

/O ^ d ) ,1 

m h^(d) a = 2 . 2A(ip.n.(-l)^ 

for n= 1,3,5, ... 

m [h^(6)-2’^. nl.i’^ erfcCg)] 

Ch_(6)-2”.nV.i’^ erfc(6)].m , 

= — .2A^p.n.(-l)^^-^ 

[hn^l(6 )-2’^'^^.(n+l) erfc(6 )] 

(aVD 2 )’^ 5 for n = 2,4,6, ... (16) 

Since typically, 6 > 1, so that 2’^.nl.i’^ erfc(6) << h^(6), the 
expression for can be written as, 

m h_(6) a - .2m App.n.(-1)’^‘^^ ( 0 :^ 02 )^ (17) 

* 

for n = 1,2,3, . . . 

Thus, 

p(t} = m C,(X (t),t) « A + E h*77X^ .2mApS.n.(-l) . 

X V H'tx 


( 18 ) 
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APPENDIX IV 


(a) Etchant for Gold : 

Potassium Iodide - 4 gm 

Iodine - 1 gm 

Deionized water — 40 c.c. 

The potassium iodide and iodine are dissolved in water 
at room temperature and stored for later use. 

(B) Etchant for Chromium : 

The etchant consists of two parts : (i) 100 gm. NaOH 
dissolved in 200 c.c. deionized water and (ii) 100 gm. 
K^FoCCN)^ dissolved in 300 c.c. deionized water. The two 
solutions are prepared and stored separately. Just before use, 
these are mixed in the ratio of 1:3 . 
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APPENDIX V 


Impurity distribution in silicon, for the case of diffusion 
in non-oxidising ambient is given by, 

/ N r , Y ax+a^D^t 

C2(x) = C^iexfci ) - e ^ .erfc( — + afo^t)] 

21/02^ 2Y" D2t 

( 1 ) 

At X - and x = X2, the eqn. (1) may be written as. 


/ V « r , V ax,+a^D,5t X 

C2(xi)= C^[erfc(— i— .)-e ^ ^ 


ayo^t 


X -1 

.erfc(--— + aY'D2t)] 


2fD^t 


( 2 ) 


ax„+a D^t 

C2(X2)= C3[erfc(3:— )-e .erfc(^^^ + o:V'D2t)](3) 


2VD,t 


In eqns. (2) and (3), and D2 are unknown. Eqns. (2) 
and (3) may be simultaneously solved for and D2 in the 
following manner. 


From eqns* (2) and (3), 


2 

X, ax«+a D«t x, C^(x, ) 

crfc(~=^ )-G ,erfc(-i + afD^t)- ' ^1 

2VD2t ^ ^2 ''^2'' 


x^ 

,[erfc(-= — ) 
2fD2t 


2 

ax,+a Dot Xrt 

- e ^ ^ .erfc(-^ + afD^t)] = 0 (4) 

2VD2t ^ 

Since C2(x3_), ^2^ ^2^* ^2 known, the eqn. (4) 

can bo numerically solved for D2 by Newton-Raphson method. With 
Dg known, the value of can be determined by substituting for 
Dg in either oqn. (2) or eqn. (3). 
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The same method as described above, has been applied to 
calculate the value of diffusion parameters in case of 

diffusion in oxidising/chloro-oxidising ambient. In this case, 

02(2) = A.erfc(~ --^^- ) (5) 

2fD2t 

In egn. (5) only A and are unknown. By considering the 
impurity concentration at two points, z^,Z 2 > the values of A and 
□2 can bo determined in a similar manner as described above. 



appendix VI 


The diffusion equation for excess self-interstitials 
(Ci-C®^) can be written as, 


D 


d( AC.) 




TT 


ACj.(eo,t) — ACj(x,o) = o 


A ^ A C. ) 


( 1 ) 


( 2 ) 


Ax ^ICl ^^^l^ 


rrAt-l /2 


*^IC1 ^l'^ 


where 


^ICl “ ^^Cl ^1 ^l^ 


Taking Laplace Transform of eqn. (1) w.r.t. 't* and writing 
AC(x,t)3 =» AC( x,s), and using the initial conditions given 
in eqn. (2), 


A Cr 

0 **** S # A C«r ^ o 


so that, 


ACj{x,s) =Aj_e'^ +A2e^ 


n k) 


(3) 


Taking Laplace Transform of the boundary conditions in eqn. (2), 


ACj(w, s) =* 0 
A tZ. 


(2a) 


*D 


i —TST * ''icr \ ‘'ici 




and 


(2b) 
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Using (2a) in (3) leads to = o. Substituting for Cj- 
in (2b) from eqn. (3), 


r r .rfS ^ n" ^ 

^ ® ^ ■♦• ^2 e I = Aj- V 


r 1 
s 


ceq 

'^ICl ” 


or 


A^ [ fsDj. + - •— £ 


Ao = [■ 


Arl/"TC 


■'loi =r 


TD. (irs+ Ifs fD..3.(Vsf 


ici, 


Vd, 


TDn 


Thus , 
AC, 


-xlTs/Tf Dj 


/> \r B— . 

Aj Y 0^ 


iTsCfst — ) 


k C®'’ 

*^101 4 

fOr 


•xlfs/VD, 


fDi 


s(^s+ :i£^) 


IfD, 


Taking inverse Laplace, 


AC, 


t’+2x ft' 

(x,t) = Aj Y ^ ^ ^ ^ Grfc(x^+ftj:^)] - 


.eq 


^ICl '"I 


fDj 

where * x/2fDj-t J 


fo, 

[ t*rfc(x )- e 

nci I Cl 


f Di t;^+2x/ 1;, 


^ ^ICl/°I 


erfc(x^+f t^)] 


For Xj^ * 0, (near the surface) 

t' 


ACUo.t) " Aj r ^ e " erfcCVt;) - [l- e " erfoClft;)] ( 4 ) 
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